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1 Introduction 



The theory of computation introduced by Blum, Shub, and Smale in [9] allows for 
computations over an arbitrary ring R. Emphasis was put, however, on the cases 
R = R OT R = C For these two cases, a major complexity result in [9] exhibited 
natural NP-complete problems, namely, the feasibility of semialgebraic or algebraic 
sets, respectively. Thus, the complexity of a basic problem in semialgebraic or alge- 
braic geometry was precisely characterized in terms of completeness in complexity 
classes. 

In contrast with discrete^ complexity theory, these first completeness results 
were not followed by an avalanche of similar results. One may say that, if NP- 
completeness exhibits a single problem with different dresses, the wardrobe of that 
problem in the real or complex settings seems to be definitely smaller than that in 
the discrete setting. 

Also in contrast with discrete complexity theory, very little emphasis was put on 
functional problems. These attracted attention at the level of analysis of particular 
algorithms, but structural properties of classes of such problems have been hardly 
studied. So far, the most systematic approach to study the complexity of certain 
functional problems within a framework of computations over the reals is Valiant's 
theory of VNP-completeness [14, 69, 72]. However, the relationship of this theory 
to the more general BSS-setting is, as of today, poorly understood. 

A recent departure from the situation above is the work focusing on complexity 
classes related with counting problems, i.e., functional problems, whose associated 
functions count the number of solutions of some decisional problem. 

In classical complexity theory, counting classes were introduced by Valiant in his 
seminal papers [70, 71]. Valiant defined as the class of functions which count 
the number of accepting paths of nondeterministic polynomial time Turing machines 
and proved that the computation of the permanent is #P-complete. This exhibited 
an unexpected difficulty for the computation of a function whose definition is only 
slightly different to that of the determinant, a well-known "easy" problem. This 
difficulty was highlighted by a result of Toda [68] proving that PH C P#P, i.e., that 
T^P has at least the power of the polynomial hierarchy. 

In the continuous setting, i.e., over the reals, counting classes were first defined 
by Meer in [50]. Here a real version #P]g of the class #P was introduced, but 
the existence of complete problems for it was not studied^. Instead, the focus of 
Meer's paper are some logical properties of this class (in terms of metafinite model 
theory). After that, in [17], an in-depth study of the properties of counting classes 
over (R, -!-,—,<) was carried out. In this setting, real computations are restricted 

^AU along this paper wc use the words discrete, classical or Boolean to emphasize that we are 
refering to the theory of complexity over a finite alphabet as exposed in, e.g., [2, 58]. 

^To distinguish between classical and, say, real complexity complexity classes, we use the sub- 
script K to indicate the latter. Also, to further emphasize this distinction, we write the former in 
sans serif. 
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to those which do not perform multiphcations and divisions. Main results in [17] 
include both structural relationships between complexity classes and completeness 
results. 

The goal of this paper is to further study #Pr (and its version over the complex 
numbers, #Pc) following the lines of [17]. A driving motivation is to capture the 
complexity (in terms of completeness results) to compute basic quantities of alge- 
braic geometry or algebraic topology in terms of complexity classes and completeness 
results. Examples for such quantities are: dimension, cardinality of 0-dimensional 
sets, geometric degree, multiplicities, number of connected or irreducible compo- 
nents, Betti numbers, rank of (sheaf) cohomology groups, Euler characteristic, etc. 
To our best knowledge, besides [17], the only known non-trivial complexity lower 
bounds for some of these quantities are in [1, 60]. For other attempts to characterize 
the intrinsic complexity of problems of algebraic geometry, especially elimination, 
we refer to [34, 48, 49]. 

Capturing the complexity of some of the above problems will help to reduce the 
contrasts wc mentioned at the beginning of this introduction. 

1. Counting classes The class #P is defined to be the class of functions 
/: {0, 1}°^ N for which there exists a polynomial time Turing machine M and a 
polynomial p with the property that for all n G N and all x G {0, 1}", f{x) counts 
the number of strings y € {0, 1}p(") such that Al accepts {x,y). 

Replacing Turing machines by BSS-machines over M in the definition above, wc 
get a class of functions / : M°° — > N U {oo}, which we denote by #P]r. Thus f{x) 
counts the number of vectors y € MP^"^^ such that M accepts {x,y). Note that this 
number may be infinite, that is, f(x) = oo. In a similar way, one defines ^P^- 

Feasibility of Boolean combinations of polynomial equalities and inequalities and 
of polynomial equations were proved to be NP]R-complete problems in [9]. These 
problems are denoted by SASr and FeaSe respectively. As one may expect, their 
counting versions ^SASr and ^FeaSr, consisting of counting the number of solu- 
tions of systems as described above, turn out to be complete in #Pr. Similarly, the 
problem #HNc consisting of counting the number of complex solutions of systems of 
polynomial equations is complete in #Pc. While we prove these results in Section 3, 
one of the goals of this paper is to show that other problems, of a basic geometric 
nature, are also complete in these counting classes. 

2. Degree, Euler characteristic and Betti numbers The study of the 
zero sets of systems of polynomial equations is the subject of algebraic geometry. 
Classically, these zero sets, called algebraic varieties, are considered in A;" for some 

algebraically closed field k. A central choice for k is k = C Civen an algebraic 
variety Z, a number of quantities are attached to it, which describe several geometric 
features of Z. Examples of such quantities are dimension and degree. Roughly 
speaking, the degree measures how twisted Z is embedded in affine space by, more 
precisely, counting how many intersection points it has with generic affine subspaces 
of a certain well-chosen dimension. Not surprisingly, an algebraic variety has degree 
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one if and only if it is an afHne subspace of C". The degree of an algebraic variety 
occurs in many results in algebraic geometry. Maybe the most celebrated of them is 
Bezout's Theorem. It also occurs in the algorithmics of algebraic geometry [24, 33] 
and in lower bounds results [16, 65]. 

The birth of algebraic topology is entangled with more than one century of 
attempts to prove a statement of Euler asserting that in a polyhedron, the number 
of vertices plus the number of faces minus the number of edges equals 2 (see [47] for 
a vivid account of this history) . A precise definition of a generalization of this sum 
is today (justly) known with the name of Euler characteristic or (justly as well) of 
Euler-Poincare characteristic. 

The Euler characteristic of X, denoted by x{^)^ is one of the most basic invari- 
ants in algebraic topology. Remarkably, it naturally occurs in many applications 
in other branches of geometry. For instance, in differential geometry, where it is 
proved that a compact, connected, differentiable manifold X has a non- vanishing 
vector field if and only if x{^) = [64, p. 201]. Also, in algebraic geometry, a 
generalization of the Euler characteristic (w.r.t. sheaf cohomology) plays a key role 
in the Riemann-Roch Theorem for non-singular projective varieties [38]. The Euler 
characteristic has also played a role in complexity lower bounds results. For this pur- 
pose, Yao [73] introduced a minor variation of the Euler characteristic. This modified 
Euler characteristic, denoted has a desirable additivity property and coincides 
with the usual Euler characteristic in many cases, e.g., for compact semialgebraic 
sets and complex algebraic varieties. 

The Euler characteristic is invariant under homotopy equivalence and the modi- 
fied Euler characteristic is invariant under homeomorphism. Thus, these quantities 
are used to prove that certain topological spaces are not homotopy equivalent or 
homeomorphic. Yet, there exist simple examples of pairs of non-equivalent spaces 
which have the same Euler characteristic. For instance the spheres and of 
dimensions 1 and 3, respectively, satisfy x{S^) = x{S^) = ^-^d they are not ho- 
motopy equivalent. A more powerful object to distinguish non-equivalent spaces is 
the sequence of Betti numbers. This is a sequence of non-negative integers hk{X), 
k > 0, associated to a topological space X, invariant under homotopy equivalence, 
and satisfying that, if the dimension of X is d, then bk{X) = for all k > d. The 
quantity bQ(X) has a very simple meaning: it is the number of connected compo- 
nents of X. Roughly speaking, for A; > 1, bk{X) counts the number of fc-dimensional 
holes of X. We have bo{S^) = bo{S^) = 1, bi{S^) = 1, bi{S^) = b2{S^) = 0, and 
^3(5'^) = 1- This shows that 5^ and 5^ are not homotopically equivalent (as one 
could well expect). The Euler characteristic and the sequence of Betti numbers are 

not unrelated. One has x(^) = 'l2k&i~^)''^li(-^) ■ 

Just as with the Euler characteristic, a version of the Betti numbers satisfying 
an additivity property was introduced by Borel and Moore [11] for locally closed 
spaces X. These Borcl-Moorc Betti numbers bf.^'\X) arc invariant under homeo- 
morphisms and arc related to the modified Euler characteristic as follows: for locally 
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closed spaces X one has x*{X) = J2k&i-'^)''b™i^)- 

3. Completeness results A semialgebraic subset of M" is defined by a Boolean 
combination of polynomial equalities and inequalities. Machines over M decide (in 
bounded time) sets which, when restricted to a fixed dimension n, are semialgebraic 
subsets of M". Therefore, this kind of sets are also the natural input of geometric 
problems in this setting. We have already remarked that deciding emptyness of a 
semialgebraic set is NPR-complete, and that counting the number of points of such 
a set is ^Pjg-complete. One of the main results in this paper is that the problem 
EuLERj^ consisting of computing the modified Euler characteristic of a semialgebraic 
set is FP]^ "-complete. The class FPj^ * is an extension of #Pr in which we allow a 
polynomial time computation with an oracle (i.e., a black box) for a function / in 
#Pig. This enhances the power of ^Pjj by allowing one to compute several values of 
/ instead of only one. 

Over the complex numbers, the situation is similar. Natural inputs for geometric 
problems are quasialgebraic sets, i.e., sets defined by a Boolean combination of 
polynomial equations. Of particular interest are algebraic varieties. We already 
remarked that deciding emptyness of an algebraic variety is NPc-complete and that 
counting the number of points of such a set is ^P^-complete. Another of the main 
results in this paper is that the problem Degree consisting of computing the degree 
of an algebraic variety is FP*^^-complete. 

The proofs of our completeness results rely on diverse tools drawn from algebraic 
geometry, algebraic topology, and complexity theory. Two of the techniques we use 
deserve, we believe, some highlight. The first one is the use of generic quantifiers, 
describing properties which hold for almost all values. A blend of reasonings in logic 
and geometry allows one to eliminate generic quantifiers in parameterized formulae. 
The basic idea behind this metliod appeared already in [36] and was used also in [7], 
but the method itself was developed in [42, 44, 45] to prove that the problem of 
computing the dimension of a semialgebraic (or complex algebraic) set is complete 
in NP]R (resp. NPc). We extend this method and use this in the completeness proofs 
of both the degree and the Euler characteristic problems. 

The second technique we want to highlight is the application of Morse theory 
for the computation of the Euler characteristic. The use of Morse functions as 
an algorithmic tool in algebraic geometry goes back to [28, 29] where the "critical 
points method" was developed to decide quantified formulae. Several algorithms to 
compute the Euler characteristic of a semialgebraic set reduce first to the case of a 
smooth hypersurface and then apply the fundamental theorem of Morse theory [3, 
13, 66]. We proceed similarly. It should be noted, however, that our reduction to 
the smooth hypersurface case is different from those in the references above since 
the latter can not be carried out within the allowed resources (polynomial time for 
real machines). 
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4. Completeness results in the Turing model In the discussion above we 
considered real solutions of systems of real polynomials and complex solutions of 
systems of complex polynomials. This coincidence between the base field for the 
space of solutions and that for the ring of polynomials used to describe solution sets 
is not necessary. While one may think of several combinations breaking it, the one 
that stands out is the consideration of real (or complex) solutions of polynomial 
systems over the integers. In practice, the difference between considering real or 
integer coefficients in the input data is reflected in the difference between the nu- 
merical analysis of polynomial systems and their symbolic computation (computer 
algebra). Note that if one restricts the input polynomials for a problem to have 
integer coefficients, then the input data for this problem can be encoded in a finite 
alphabet and may be considered in the classical setting. To distinguish this dis- 
cretized version from its continuous counterpart we will add a superscript "0" in the 
problem's name. Thus, for instance, HNc'^ is the problem of deciding the existence 
of complex solutions of a system of integer polynomial equations and ^HNc^ is the 
problem of counting the number of these solutions. 

The complexity of computer algebra algorithms for, say, HNc° is described using 
discrete models of computation (e.g., Turing machines). For instance, relatively 
recent results [24] show that HNc*^ E PSPACE, and an even more recent result 
of Koiran [40] shows that, assuming the generalized Riemann hypothesis, HNc'^ G 
RP'^'^. On the other hand, it is well-known (and rather trivial) that HNc° is NP- 
hard. The complexity of problems like Fea^ or SAS^ is much less understood, the 
gap between their known lower NP and upper PSPACE bounds being much larger. 

In this paper we introduce two new counting complexity classes in the discrete 
setting namely, GCC and GCR. These classes are closed under parsimonious reduc- 
tions and located between #P and FPSPACE. The problem tt^HNc'^ is complete in 
GCC and the problems t^SASk" and ^^FeaSk^ are complete in GCR. In addition, 
we also prove that Degree" and Euler^° are complete in FP^^*- and FP^^'^, re- 
spectively, and that EULER^, the problem of computing the (non- modified) Euler 
characteristic of a basic semialgebraic set, is complete in FP*^^'^. 

Canny [18] showed that the problem #CCg of counting the number of connected 
components of a semialgebraic set described by integer polynomials is in FPSPACE. 
On the other hand, a result by Reif [60, 61] stating the PSPACE-hardness of a 
generalized movers problem in robotics easily implies the FPSPACE-hardness of the 
problem #CC^ 

We give an alternative proof of the FPSPACE-hardness of #CC^ following the 
lines of [17]. Extending this, we prove that the problem BETTl(fc)^ of computing 
the kth Betti number of the real zero set of a given integer polynomial is FPSPACE- 
hard, for fixed k £ N. We also prove that the problem BM-Betti(A:){^ of computing 
the fcth Borel-Moore Betti number of the set of real zeros of a given integer poly- 
nomial is FPSPACE-hard. Note that, for k > 1, the membership of Betti(A:)^ and 
BM-Betti(A;)^ to FPSPACE is, as of today, an open problem. 
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State-of-the-art algorithmics for computing the Euler characteristic or the num- 
ber of connected components of a semialgebraic set suggests that the former is 
simpler than the latter [3, 4]. In a recently published book [5, page 547] it is ex- 
plicitly observed that the Euler characteristic of real algebraic sets (which is the 
alternating sum of the Betti numbers) can be currently more efficiently computed 
than any of the individual Betti numbers. 

Our results give some explanation for the observed higher complexity required for 
the computation of the number of connected components (or higher Betti numbers) 
compared to the computation of the Euler characteristic. Namely, EULER^ is FP^^*^- 
complete, while BETTi(fe)^ is FPSPACE-hard. Thus the problem Betti(/c)^ is not 
polynomial time equivalent to EULER^ unless there is the collapse of complexity 
classes FP^^*^ = FPSPACE. 

A similar observation for the Euler characteristic and the Betti numbers in the 
context of semi-linear sets and additive machines was made in [17, Corollary 5.23]. 

5. Organization of the paper We start in Section 2 by recalling basic facts 

about machines and complexity classes over M and C as well as about semialgebraic 
and algebraic sets. Then we define in Section 3 the counting complexity classes 
#Pk and #Pc, introduce different notions of reduction, and prove some basic com- 
pleteness results. The technique of generic quantifiers is described in Section 4 and 
then used in Section 5 to prove the completeness result for DEGREE. The proof 
of this result is preceded by the exposition of some basic facts about smoothness 
and transversality, which lead to a concise way of expressing the degree by a pa- 
rameterized first order formula. We prove the completeness of EULER^ in Section 7 
after recalling some basic facts from algebraic and differential topology in Section 6. 
Section 8 deals with complexity in the discrete setting. We define the classes GCC 
and GCR and, besides some basic completeness results, we prove the completeness 
of Degree'^ in GCC and of Euler^ and Euler^° in GCR. Finally, we prove the 
FPSPACE-hardness of the problems Betti(A;)^ and BM-Betti(A;)}^. Wc close the 
paper in Section 9 with a summary of problems and results, and with some selected 
open problems in Section 10. 

Acknowledgment. We are thankful to Saugata Basu and Pascal Koiran for 
helpful discussions while writing this paper. We are specially indebted to Pascal 
Koiran since many of his results have been a great source of inspiration for us. 

2 Preliminaries about real machines 

2.1 Machines and complexity classes 

We denote by M°° the disjoint union R°° = |Jn>o^"' where for n > 0, is the 
standard n-dimensional space over M. The space is a natural one to represent 
problem instances of arbitrarily high dimension. For x G C M°°, we call n the 
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size of X and we denote it by size(a;). Contained in M°° is the set of bitstrings {0, 1}°° 
defined as the union of the sets {0, l}", for n G N. 

In this paper we wih consider BSS-machines over M as they are defined in [8, 9]. 
Roughly speaking, such a machine takes an input from M°°, performs a number 
of arithmetic operations and comparisons following a finite list of instructions, and 
halts returning an element in ]R°° (or loops forever). 

For a given machine M, the function (pM associating its output to a given input 
X G is called the input-output function. We shall say that a function / : M°° — > 
MJ^, A; < oo, is computable when there is a machine M such that / = tpM- Also, a set 
A C is decided by a machine M if its characteristic function XA '■ — ^ {Oj 1} 
coincides with ipM- So, for decision problems we consider machines whose output 
space is {0, 1} C M. 

We next introduce some central complexity classes. 

Definition 2.1 A machine M over R is said to work in polynomAal time when there 
is a constant c G N such that for every input x G M°^, M reaches its output node 
after at most size(a:;)'^ steps. The class Pr is then defined as the set of all subsets of 
M°° that can be accepted by a machine working in polynomial time, and the class 
FP]R as the set of functions which can be computed in polynomial time. 

Definition 2.2 A set A belongs to NPu if there is a machine M satisfying the 

following condition: for all x G x G A iff there exists y G M°° such that M 
accepts the input {x, y) within time polynomial in size(a;). In this case, the element y 
is called a witness for x. 

Remark 2.3 (i) In this model, the element y can be seen as the sequence of guesses 
used in the Turing machine model. However, we note that in this definition no 
nondeterministic machine is introduced as a computational model, and non- 
determinism appears here as a new acceptance definition for the deterministic 
machine. Also, we note that the length of y can be easily bounded by the time 
bound p(size(x)). 

(ii) Machines over C arc defined as those over M. Note, though, that branchings 
over C are done on tests of the form zq = 0. The classes Pc, NPc, etc., are 
then naturally defined. 

In [8, Chapter 18] models for parallel computation over M are defined. Using 
these models, one defines PARjr to be the class of subsets of M°°, whose characteristic 
function can be computed in parallel polynomial time. Also, one defines FPAR]r to 
be the class of functions computable in parallel polynomial time such that size(/(x)) 
is bounded by a polynomial in size(x). 
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2.2 Algebraic and semialgebraic sets 

Algebraic geometry is the study of zero sets of polynomials (or of objects which 
locally resemble these sets). Standard textbooks on algebraic geometry are [31, 54, 
63]. For information about real algebraic geometry we refer to [6, 10]. 

We very briefly recall some definitions and facts from algebraic geometry, which 
will be needed later on. 

An algebraic set (or affine algebraic variety) Z is defined as the zero set 

Z = Zih, ...,/,):= {x G C" I /i(x) = 0, ... , fr{x) = 0} 

of finitely many polynomials /i, . . . , € C[Xi, . . . , Xn]- The vanishing ideal 2{Z) 
of Z consists of all the polynomials vanishing on Z. Note that I{Z) might be strictly 
larger than the ideal / generated hy fi, . . . , fr- Actually, by Hilbert's Nullstellensatz, 
can be characterized as the so-called radical of the ideal /. 

A usual compactification of the space C" consists of embedding C" into P"'(C), 
the projective space of dimension n over C. Recall, this is the set of complex lines 
through the origin in C""^^ and ^ P"(C) maps a point a; G to the line 
in C"+^ passing through the origin and through (l,x). The notion of an affine 
algebraic variety extends to that of a projective variety by replacing polynomials 
by homogeneous polynomials in C[Xo, Xi, . . . , Xn], for which elements of P"'(C) are 
natural zeros. The embedding C" ^ P'*(C) extends to the algebraic subsets of C" 
by defining, for any such set Z, its projective closure Z as the smallest projective 
variety in P"(C) containing Z. 

A basic semialgebraic set S C R"' is defined to be a set of the form 

S = {x G M'^ I g{x) = 0,h{x) > 0, . . .,fr{x) > 0}, 

where g, fi, . . . , fr are polynomials in M[Xi, . . . , X„]. We say that iS C is a 
semialgebraic set when it is a Boolean combination of basic semialgebraic sets in 
M". Every semialgebraic set S can be represented as a finite union S = SiU . . .U St 
of basic semialgebraic sets.^ 

We will consider algebraic or semialgebraic sets as input data for machines over M 
or C. These sets are encoded by a family of polynomials describing the set as above. 
To fix ideas we will assume, unless otherwise specified, that semialgebraic sets are 
given as unions of basic semialgebraic sets. So, properly speaking, the input data is 
not the set itself but a description of it. Also, we have to define how polynomials 
themselves are encoded as vectors of real (or complex) numbers. However, it will 
turn out that our results have little dependence on the choice of the representation 
of the semialgebraic set and on the encoding of the polynomials, cf. Remark 9.1. 

A polynomial / = X^gg/ Ue x^^ • • • x^" is represented in the sparse encoding by a 
list of the pairs {ug, e) for e G i", where / = {e G | -Ug 7^ 0}. The coefficients 

^This respresentation is said to be in Disjunctive Normal Form. A representation in Conjuntive 
Normal Form is defined in the obvious manner. 
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Ue are given as real (or complex) numbers, while the exponent vector e is thought 
to be given by a bit vector of length at most 0(nlogdcg/). Let |/| be the total 
number of terms and S := max{2,deg/}. Then size(/) := \I\nlog5 is defined to 
be the sparse size of /. The sparse size of a set of polynomials /i, . . . , is defined 
as size(/j). To fix ideas, we will always assume that polynomials are given 

by the sparse encoding. If we are dealing with integer polynomials /, we will also 
consider their sparse bit size, which is defined as the sparse size of / multiplied by 
the maximum bit size of the occuring integer coefficients. 

We remark that another way of encoding polynomials is the dense encoding. 
Here, a polynomial of degree d in n variables is given by the list of its {''Y) coefH- 
cients, which has therefore the size of this combinatorial number. Yet another way 
is to encode the polynomial by a straight-line program computing it, cf. [8, 16]. In 
this case, the size of the encoding of / is the length of the straight-line program. 

2.3 Some known completeness results 

We first recall the basic notions of reduction for classes of decision problems.^ 

Definition 2.4 1. Let S,T C R°°. We say that M°° ^ is a reduction 
from S to T if can be computed in polynomial time and, for all x G M°°, 
X G -S if and only if ip{x) G T. 

2. We say that S Turing reduces to T if there exists an oracle machine which, 
with oracle T, decides S in polynomial time. 

3. Let C be any class of subsets of R°°. We say that a set T is hard for C if, for 
every S ^ C, there is a reduction from S to T. We say that T is C-complete if, 
in addition, T & C. 

4. The notions of Turing-hardness or Turing-completeness are defined similarly. 

The extension of this definition to C is immediate. 

The following problems describing variants of the basic feasibility problem over 
M and C were introduced and studied in [9]. 

HNc {Hubert's Nullstellensatz) Given a finite set of complex multivariate poly- 
nomials, decide whether these polynomials have a common complex zero. 

FeaSr {Polynomial feasibility) Given a real multivariate polynomial, decide whether 
it has a real root. 

SAS]K {Semialgebraic satisfiability) Given a semialgebraic set 5, decide whether 
it is nonempty. 

''This definition is actually for a class C containing NPe fl coNPg. To define PM-completeness, a 
stronger notion of reduction is necessary. 
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In [9], the following fundamental completeness result was proved. 

Theorem 2.5 The problem HNc is NPc-compJete and the problems FeaSr and 
SASr are NPu-compJete. □ 

Consider the following decision problems related to the computation of the di- 
mension of algebraic or semialgebraic sets. 

DiMc {Algebraic dimension) Given a finite set of complex polynomials with affine 
zero set Z and d G N, decide whether dimZ > d. 

DiMr (Semialgebraic dimension) Given a semialgebraic set S and d G N, decide 
whether dim 5 > d. 

We denote by DlM^ the restriction of the problem DiMc to input polynomials 
with integer coefficients. This problem can be encoded in a finite alphabet and may 
thus be studied in the classical Turing setting. The problems DiM^ and HN£ are 
defined similarly. 

Koiran [42, 45] significantly extended the list of known geometric NPc- or NFr- 
complete problems by showing the following. 

Theorem 2.6 (i) DiMc is NPc-compiete, and DlM^. is equivalent to HNc° with 
respect to polynomial-time many-one reductions. 

(ii) DiMjR is NPiR-compiete, and DlM^ is equivalent to FeaSir° with respect to 
polynomial-time many-one reductions. □ 

3 Counting Complexity Classes 

Definition 3.1 We say that a function /: M°° — > NU{oo} belongs to the class #P]k 
if there exists a polynomial time machine M over R and a polynomial p such that, 
for all a; G M*^, 

f{x) = \{y€ M^'(") I M accepts {x,y)}\. 

The complexity class FP£" * consists of all functions /: which can be 

computed in polynomial time using oracle calls to functions in 

Remark 3.2 (i) The class #Pig is the one defined by Meer in [50]. 

(ii) The counting classes #Pc and FP*^^ are defined mutatis mutandis. Also, 
replacing M by Z2 in Definition 3.1 one obtains the classical #P. 

We next define appropriate notions of reduction and completeness. 
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Definition 3.3 1. Let f,g: R°° ^ N Li {oo}. We say that R°° ^ R°° is a 
parsimonious reduction from f to g if tp can be computed in polynomial time 
and, for all x G R°°, f{x) = g{(fi{x)). 

2. We say that / Turing reduces to g if there exists an oracle machine which, 
with oracle 5, computes / in polynomial time. 

3. Let C be #P]r or FP^ We say that a function g is hard for C if, for every 
f € C, there is a parsimonious reduction from / to g. We say that g is 
C-complete if, in addition, g G C. 

4. The notions of Turing-hardness or Turing-completeness are defined similarly. 
The extension of this definition to C is immediate. 

We define now the following counting versions of the basic feasibility problems 
HNcFeaSr, and SASr. 

T^HNc {Algebraic point counting) Given a finite set of complex multivariate poly- 
nomials, count the number of complex common zeros, returning 00 if this 
number is not finite. 

^FeaSr {Real algebraic point counting) Given a real multivariate polynomial, 
count the number of its real roots, returning 00 if this number is not finite. 

#SAS]K {Semialgebraic point counting) Given a semialgebraic set S, compute its 
cardinality if S is finite, and return 00 otherwise. 

As was to be expected, these counting problems turn out to be complete in the 
classes and #Pr, respectively. In the sequel, given n G N, we denote by [n] the 
set {!,..., n}. 

Theorem 3.4 (i) The problem #HNc is ^F^-complete (with respect to parsimo- 
nious reductions). 

(ii) The problems #FeaS]k and #SAS]k are i^P^-complete with respect to Taring 
reductions. 

Proof. For part (i) simply check that the reductions given in the corresponding 
NPc-completeness result by Blum, Shub and Smale [9] (see also [8]) are parsimonious. 

The proof of part (ii) requires a more careful look at the reduction in [8] . In this 
proof, a machine M solving a given problem in NPr is considered and a reduction is 
established, which associates to every input u> G R°°, a conjunctive normal form ijj^ 
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(the fact that there is only one equahty in each clause is achieved by adding squares). 
An important point to remark here is that, while the cardinality of / is bounded by 
a polynomial in the size of a;, the cardinalities of the sets Jj are independent of u 
and depend only on M. 

Now consider one of the clauses of 

gi{x) = y V fijix)>0. (1) 

Considering that gi may be at a point x either = or 7^ 0, and that fij may be 
either < 0, = or > we have 2 x 3^' possibilities for the signs of gi, fa, ■ ■ ■ fin 
a point X. Prom them, only = 2 x 3''* — 2^' satisfy the clause (1). We conclude 
that we can rewrite this clause as an exclusive disjunction of Ki conjunctions of the 
form 

gi{x)/\iQ A l\ fij{x)UijQ (2) 

where Aj G {=, 7^} and G {<, =, >}. Now replace in (2) the occurrences 

gi{x) / by gi{x)zi -1 = 0, 

/ij(x)>0 by /ij(x)|/| - 1 = 0, 

hj{x)<Q by /,^.(x)4 + i = o, 

fi, (x) = by (x) = A 4 -1 = 0. 

This yields a system of equalities which has, for every solution x of (2), exactly 
2''» solutions in the variables x,y,z. Now, for £ G [Ki], reduce the system in (2) 
corresponding to ^ to a single equation Fii{x, y,z) =0 by adding squares and the 
clause (1) to an equation F*{x,y,z) = by taking F* = Y[f=i^ii- Note that, 
for each solution x of there are exactly 2*" different solutions {x,y,z) of the 
polynomial 

F:=F^{x,y,zf + --- + F:;,{x,y,zf 



where Tn is the cardinality of / and v — t± -\- . . . -\- Vj-, 

The #P]K-Turing-hardness of FeaSk now follows. Finish the reduction above by 
quering FeaSk for the polynomial F and divide the result by 2^^. □ 



Remark 3.5 The proof of Theorem 3.4 shows that the version of SASr with semi- 
algebraic sets given in conjunctive normal form is #PR-complete with respect to 
parsimonious reductions. 

Proposition 3.6 If f € #Pr then, for all x € for which /(x) is finite, the 
bit-size of /(x) is bounded by a polynomial in the size of x. 
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Proof. To prove the statement note that, given x G there exist polynomials 
q such that the set of witnesses for x is a semialgebraic subset of M^^") defined by 
a union of at most 2'^'^"'' basic semialgebraic sets, each of them described by a system 
of at most g(n) inequalities of polynomials in p{n) variables with degree at most 
2g(n) jf ^]jig gg^ finite, its cardinality coincides with the number of its connected 
components. Now use the bounds for the number of connected components of such 
basic semialgebraic sets (see e.g. [16, Thm. 11.1] or [8, Prop. 7, Chapt. 16]), which 
follow from the well-known Olemik-Petrovski-Milnor-Thom bounds [52, 56, 57, 67]. 

□ 

We next locate the newly defined counting complexity classes within the land- 
scape of known complexity classes. 

Theorem 3.7 We have FPj^ * C FPARk. (To interpret this, represent oo by an 
element of M - N.) 

Proof. By Theorem 3.4(i), it is sufficient to prove that #SASk belongs to FPARjr. 
By Theorem 2.6(ii), the problem of computing the dimension of a semialgebraic set 
IS m FP^^, and therefore, in FPARr. We use this to compute #SASk as follows. 
Given a semialgebraic set, we check whether it is zero dimensional. If yes, we return 
its number of connected components, otherwise we return oo. This is in FPARr due 
to the main result in [4, 30, 35]. □ 

Remark 3.8 Versions of Proposition 3.6 and of Theorem 3.7 hold over C as well, 
with proofs similar to those over R. 

The following lemma will be useful later on. It is an immediate consequence of 
the definition of the counting classes. 

Lemma 3.9 Let f: R°° x {0, 1}°° -> N be a function in #Fj^. Assign to f and a 
polynomial p the following function g: R°° — > N obtained by summation: forx G M", 

2/e{o,i}J'(") 

Then g belongs to #Pir. A similar statement holds over C. □ 

4 Generic quantifiers 

Our completeness results for Degree and EULER^ crucially depend on Koiran's 
method [42, 44, 45] to eliminate generic quantifiers in parameterized formulas. In 
this section, we further develop Koiran's method in order to adapt it to our purposes. 

The main difference to [42, 44, 45] is the introduction of the notion of a partial 
witness sequence (compared to the notion of a witness sequence from [42]). 
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4.1 Efficient quantifier eUmination over the reals 

For convenience of the reader, we recall a well-known result about efficient quantifier 
elimination over the reals from Renegar [62, Part III]. In the sequel denotes the 

set of first order formulas over the language of the theory of ordered fields with 
constant symbols for real numbers. The subset of formulas with constant symbols 
for and 1 only, is denoted by 

Theorem 4.1 Let F be a formula in in prcncx form with k free variables, 
n bounded variables, w alternating quantifier blocks, and m atomic predicates given 
by polynomials of degree at most S >2 with integer coefficients of bit-size at most i. 
That is, F has the form 

(Qix(i) eM"i)...(Q^x('^) GM"-)G(y,xW,...,x(^)) 

with alternating quantifiers Qi € {3, V} and free variables y = {yi, . . . , y^) G M*^; the 
quantifier free formula G is a Boolean function of m atomic predicates 

5,(y,x«,...,x("'))A,0, l<j<m, 

where the gj are integer polynomials of degree at most 6 and with coefficients of 
bit-size at most i. Hereby, Aj is any of the standard relations {>,>,=, 7^, <,<}. 
Tiien F is equivalent to a quantifier-free formula F' in disjunctive normal form 

I Ji 

1=1 j=i 

where hij are integer polynomials with degree at most D and bit-size at most L, 
and such that 

w w 

log D < 2°("') Yl ^i) ^og{mS), log L < 2°^'^^ JJ \og{m5) + log(A; + 1). 

i=l 1=1 
Moreover, the number M := Yli=i sitomic predicates satisfies the bound 

w 

logM < 2^("')fc(^ JJni) \og{mS). 

i=l 

4.2 Construction of generic points 

Definition 4.2 Let F G have free variables ai, . . . , Ofc. We say that F is Zariski- 
generically true if the set of values a G not satisfying F{a) has dimension strictly 
less than k. We express this fact by writing V*aF(a) using the generic universal 
quantifier V*. 
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Remark 4.3 (i) Let F G have k free variables and coefficient field K, i.e., K 
is the field generated by the coefficients of all the polynomials occuring in F. 
Then \/*aF{a) is equivalent to each of the following statements: 

(a) {a G R*^ I F{a)} is dense in M.'^ with respect to the Euclidean topology, 

(b) Ve G R Va e M'^ 3a' eW'' (e > ^ F{a') A \\a - a'\\ < e), 

(c) Va G M'^( ai, . . . ,afe algebraically independent over K F{a)). 

Part (b) shows that V*a -F(a) can be expressed by a first order formula. Hence 
by using the generic quantifier we still describe semialgebraic sets. 

(ii) One can define the generic existential quantifier 3* by 

3*aF{a) = ^V*a^F(a). 

Note that 3*aF{a) iff the set of values a G satisfying F{a) has dimension k. 
We may say that F is Euclidean- generically true. 

(iii) For first order formulas over the language of the theory of fields with 
constant symbols for complex numbers, one can define V* and 3* just as above. 
It is not difficult to see, however, that these two quantifiers coincide over C. 
That is, Zariski genericity equals Euclidean genericity. 

The following result was proved in Koiran [45, Cor. 1]. 

Proposition 4.4 (i) Let F G #"1^ be in prenex form with k free variables, n 
bounded variables, w alternating quantifier blocks, and m atomic predicates 
given by polynomials of degree at most S >2 with integer coefficients of bit-size 
at most I. If F is Zariski-generically true, then a point a G Z'^ satisfying F can 
be computed by a division-free arithmetic straight-line program T of length 
0{kn'^ \og{m5) + log£) having 1 as its only constant and no inputs. 

(ii) There exists a Turing machine which, with input {k,n,w,m,6,l), computes T 
in time polynomial in the length ofT. This machine docs not depend on F. 

Since we will need the proof method behind this result later on, we recall the 
proof. A first ingredient is the following easy lemma, whose proof can be found for 
instance in [41]. 

Lemma 4.5 For positive integers k, L, D recursively define 

ai := 2^, aj := 1 + ai(D + l^^^af.i for 2<j<k. 

Then h{ai, . . . , a^) / for any integer polynomial h in k variables of degree at 
most D and coefficients of absolute value less than 2^. □ 
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The sequence ai, . . . , in Lemma 4.5 can be computed by a straight-line pro- 
gram r performing 0(A;log-D -|- logL) arithmetic operations and which has 1 as its 
only constant. 

Proof of Proposition 4-4- Put S" = {a € | F{a) holds}. We use Theo- 
rem 4.1 to replace the formula F by an equivalent quantifier free formula F' = 
Vi=i A/=i hijAijO and claim that 

P|{aGM'= I hij{a)^0} C S. 

Otherwise, there would be some a € M'^' — such that hij{a) ^ for all Since 
the sign of hij does not change in a small neighborhood of a, M*^ — 5 would contain 
some ball around a. And this contradicts the assumption that S is dense in W^. 

Let D and L be the upper bounds on the degree and bit-size of the polynomials 
occurring in F' , given by Theorem 4.1. According to Lemma 4.5, we can compute 
a point a ^ iJ^ satisfying hij{a) / 0, for all i,j and thus F{a), by a straight-line 
program with 0{k\ogD + logL) arithmetic operations. By plugging in the bounds 
on D and L the claim follows (use (rij^i)''^^"' ^ n/w). □ 

4.3 Partial witness sequences 

Let ii' C M and a G M*^ with components algebraically independent over K. By 

Remark 4.3(i)(c), for any formula F with coefficient field contained in the im- 
plication (V*aF(a)) F{a) holds. Thus a may be interpreted as a partial witness 
for V*aF(a). 

Remark 4.6 The converse of the implication above, i.e., F{a) =^ (V*aF(a)) does 
not hold in general. Actually, a point a G M'^ as above and such that F[a) is true 
only ensures Euclidean genericity: we have F{a) (3*aF(a)). 

Over C, the equivalence (y*aF{a)) 4^ F{a) holds since Euclidean and Zariski 
genericity are equivalent. 

Given a formula F{u, a) we are now interested in partial witnesses for its Zariski- 
genericity property which can be used for all values of the parameter u. This may 
not be attainable with a single partial witness, but it turns out to be doable by 
using short sequences of such witnesses and taking a majority vote. Recall that [n] 
denotes the set {1, . . . , n}. 

Definition 4.7 Let F{u,a) G J^k with free variables u ^ W and a G M'^. A 
sequence o: = (cci, . . . ,Q:2p+i) G {M}^)^'^^'^^^ is called a partial witness sequence for 
F if 

Vn G W {y*a G F{u,a) =^ \{i G [2p + 1] \ F{u,ai)}\ > j)). (3) 
We denote the set of partial witnesses of F by PW{F). 
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Lemma 4.8 PW{F) is Zariski dense in R'=(2p+i). 



Proof. The proof is by a transcendence degree argument similar as in [42, Thm. 5.1]. 
Let K be the coefficient field of F. We interpret (3) as a first order formula in 
with free variables oti,. . . , a2p+i and coefficient field K. Applying Remark 4.3(i)(c) 
to this formula, it is enough to show that a G PW{F) for any cx e Rfe(2p+i) with 
components algebraically independent over K. 

Take such a and let u eW. Let K' be the field extension of K generated by the 
components of u and let K" be the field extension of K' generated by the components 
of a. Then the transcendence degree of K" over K' is at least k{2p + 1) — p. Let B 
be a transcendence basis of K" over K' consisting of components of a. Then B can 
omit components of at most p of the Qj's. The remaining Qj's have algebraically 
independent components over K' and therefore F{u, ctj) holds true for them. Thus 
\{ie[2p + l]\F{u,ai)}\>p. □ 

The next theorem is similar to [45, Thm. 3]. 

Theorem 4.9 (i) Let F{u, a) G be in prenex form with free variables u 

and a E M.'^, n bounded variables, w alternating quantiser blocks, and m 
atomic predicates given by polynomials of degree at most S > 2 with integer co- 
efficients of bit-size at most £. Thenapointo. G PVF(F)nZ'^'(^P+^) can becom- 
puted by a division-free straight-line program T of length (kp)'-^^^^ n'" \og{m6) + 
O(log^) having 1 as its only constant and no inputs. 

(ii) There exists a Turing machine which, with input {p, k, n, w, m, S, i), computes T 
in time polynomial in the length ofF. This machine does not depend on F. 

Proof. We first replace the formula F by a quantifier free formula F' according to 
Theorem 4.1. Let M be the number of atomic predicates of F', and D and L be the 
degree and the bit-size of the occuring polynomials, respectively. We have 

logD < 0(n'"log(m5)), logL < Oin"" logimS) +log{p + k + £)), 

and 

logM < 0{kn'"log{m6)). 

We replace the generic quantifier in formula (3) according to Remark 4.3(i)(b) and 
thus write the formula as 

W Ve Va 3a' < V {F'{u, a') A \\a - a'\\ <e) /\ F'{u, a^)) , 

where / runs over all p+ 1-element subsets of [2p+ 1]. This formula, let us call it ip, 
defines PW{F) and is therefore Zariski-generically true by Lemma 4.8. We may 
therefore apply Proposition 4.4 to the prenex formula ijj. Note that ^ has k{2p-\- 1) 
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free variables and 2k+p+l bounded variables, two quantifier blocks, and polynomials 
of degree at most D and bit-size at most L. The number of atomic predicates of tp 
equals (2p + 2)M + 2. Proposition 4.4 therefore implies that we may compute an 
integer point in PW{F) by a straight-line program with 0{kp{k -|- p)^ log(MZ)) -|- 
logL) arithmetic operations. The latter can be bounded by (/cp)'-'^^) n'"log(m5) -|- 
O(log^). This shows part (i). Part (ii) follows from part (ii) in Proposition 4.4. □ 

Remark 4.10 (i) It follows from part (ii) Theorem 4.9 that the element a in 
part (i) of this theorem can be computed by a machine over M or C, upon 
input (p, fc, n, u;, m, (5, t), in time order of the length of V. Note, however, that 
this computation may not be possible within these time bounds in the classical 
setting since the bit-size of the components in a. grows exponentially fast. 

(ii) Over the field C one can define the stronger notion of witness sequence. For this 
we replace in formula (3) of Definition 4.7 the implication from left to right 
by an equivalence. The analogue of Lemma 4.8 is then true and therefore 
witness sequences can be computed by "short" straight-line programs as in 
Theorem 4.9. This approach was taken in Koiran [42] to devise a method to 
compute dimensions of algebraic sets in NPc. Over the reals one cannot work 
with witness sequences, but the method can be saved by working with partial 
witness sequences as described above. 

5 Complexity of the geometric degree 

The (geometric) degree degZ of an algebraic variety Z embedded in affine or pro- 
jective space can be interpreted as a measure for the degree of nonlinearity of Z. 
A detailed treatment of this notion can be found in standard textbooks on alge- 
braic geometry [31, 54, 63]. In this section "dimension" always refers to complex 
dimension. 

Definition 5.1 Let Z C C"^ be an algebraic set of dimension d > 0. If Z is 
irreducible then its (geometric) degree deg Z is the number of intersection points of 
Z with a generic affine subspace of codimension d. If Z is reducible then its degree 
is the sum of the degrees of all irreducible components of Z of maximal dimension.^ 
The degree of the empty set is defined as 0. 

We are going to study the following problem in the computational model of 
machines over C. 

Degree ( Geometric degree) Given a finite set of complex multivariate polynomi- 
als, compute the geometric degree of its affine zero set. 

^We note here that in algebraic complexity it is common to define the degree of a reducible 
variety as the sum of the degrees of all irreducible components (cf. [16]). 
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Here is the main result of this section. 

Theorem 5.2 The problem DEGREE is FPj^ ^ -complete for Turing reductions. 

The difficult ^art of the proof is the upper bound, i.e., the membership of 
Degree to FP^ To show this membership, we have to describe a polynomial 
time algorithm over C, which computes the degree using oracle calls to #Pc. The 
basic idea of our Degree algorithm is very simple. Let /i, . . . , be an instance for 
Degree and denote its zero set by Z. We first compute the dimension d = dimZ 
by calls to HNc-oracles using Theorem 2.6. By definition, degZ is the number of 
intersection points of Z with a generic affine subspace A of codimension d. If we 
could compute such an A, then the number of intersection points could be obtained 
by a call to #HNc. 

The difficulty is how to compute a generic affine subspace. Of course, the obvious 
way to turn this idea into an algorithm would be to choose the subspace A at random. 
This would yield a randomized algorithm for computing the degree. However, our 
goal is to choose A deterministically. We will do so using partial witness sequences 
for parametrized formulas as described in Section 4, for which we need to concisely 
express the degree. If Aa denotes an affine subspace of C" of codimension d encoded 
by the parameter a e C^, then we have by the definition of degree 

V*aGC'^ \ZnAa\ =degZ. (4) 

It is clear that the above statement can be expressed by a first-order formula over C. 
However, the obvious way to do this leads to a formula with exponentially many 
variables since deg Z can be exponentially large. 

Our goal is thus to express (4) in a more concise way. This will be achieved by 
using the notion of transversality (see Lemma 5.6). However, the translation of the 
transversality condition into a concise first order formula is a little subtle and will 
require some further ideas (see Lemma 5.9). 

5.1 Smoothness and transversality 

An important notion in algebraic geometry is that of a smooth point in a variety. 
To define smoothness we use Zariski tangent spaces. 

Definition 5.3 Let Z C C" be an algebraic set, x E Z, and /i, . . . , be generators 
of the vanishing ideal I{Z) of Z. The Zariski tangent space T^Z of Z at a; is defined 

by 

TxZ = Z{d^fi, . . . ,dxfr) 

where the differential of f at x, dxf'-C" — > C, is the linear function defined by 
dxfX = Yl^=i 9xjf{x)Xj. We say that a; is a smooth point of Z if the dimension of 
TxZ equals the local dimension dim^; Z o{ Z &t x. A point in Z which is not smooth 
is said to be a singular point of Z. 
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Remark 5.4 Note that T^Z is easy to compute from a set of generators of 2r(^), 
but it may not be so, if instead we only have at hand an arbitrary set of polynomials 
with zero set Z. 

Definition 5.5 Let 2' C C" be an algebraic set of dimension d and ^ C C" be an 
affine subspace of codimension d. 

1. ^ is called transversal to Z at x £ Z D A iS x is a smooth point of Z and 
T^Z®T^A = C". 

2. We say that A is transversal to Z when A is transversal to Z at all intersection 

points X € Z n A and if, additionally, there are no intersection points of Z 
and A at infinity. No intersection points at infinity means that Z Ci A C. C", 
where Z and A are the projective closures in P"(C) of Z and A. 

In the following, wc will parametrize affinc subspaccs of codimension d, as follows. 
We denote by Aa C C" the affine subspace of C" described by the system of linear 
equations gi{x) = 0, . . . ,gd{x) = with coefficient vector a G C'*, where h = 
d{n + 1) = 0{'n?'). Note that dim^l > n — d for all a and V*a dim = n — d. 

The following lemma shows that the transvcrsality of A to Z can be used to 
certify that the number of intersection points of Z and A equals degZ. 

Lemma 5.6 If Z C C" is an algebraic set of dimension d and h = d{n + 1), then 
we have: 

(i) V*a G C'^ Aa is transversal to Z 

(ii) Va G C'' {Aa is transversal to Z ^ \Z f] Aa\ = degZ). 

Proof. This lemma is proved in Mumford [54, §5 A] for irreducible projective vari- 
eties Z. It remains to show that it extends to the case where Z is affine and reducible. 
Let Zi, . . . ,Zt be the irreducible components of Z. A dimension argument shows 
that for a generic a, Aa does neither meet the components Zi of dimension less 
than d, nor the intersections Zi n Zj for i < j. Similarly, Aa does not meet Zi — Zi 
for generic a. Hence (i) follows from the corresponding statement for irreducible 
projective varieties. 

For proving (ii) we assume that Aa is transversal to Z . Then codXvuAa = d and 
each point x G Z fl Aa is a smooth point of Z of local dimension d. Hence there is 
exactly one irreducible component of Z passing through x and this component has 
dimension d. We therefore have \Z f] Aa\ = J2i=i \^i^-^a\ where Zi, . . . , denote 
the irreducible components of dimension d. Moreover, Aa is transversal to each of 
these Zi, hence \Zi Aa\ = degZj by [54, §5A, Thm. 5.1]. Altogether, we obtain 
\ZnAa\ = J2i=i deg Zi = deg Z by the definition of the degree of reducible algebraic 
sets. □ 
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5.2 Expressing smoothness and transversality 

Lemma 5.6 suggests to use transversality to concisely express degree. But, in turn, 
to express transversality a difficulty may arise. When we try to describe the Zariski 
tangent space of Z at a point x, the given equations /i = 0, . . . , /r = for Z might 
not generate the vanishing ideal of Z, since multiplicities might occur. In other 
words, the ideal generated by fi,. . . ,fr might be different from the radical ideal, 
and it is not clear how to compute generators of the radical within the resources 
allowed. As a way out, we will express the tangent space and the transversality 
condition at x by a first order formula, in which all information regarding Z is given 
by a unary predicate expressing membership of points to Z.^ 

To do so we will use the notion of intersection multiplicity, so we next recall 
some facts about it. For more on this, the book by Mumford [54] is an excellent 
reference fitting well our geometric viewpoint.^ 

Definition 5.7 Assume that Z C C"^ is an irreducible variety of dimension d and 
let Aa C" be an affine subspace of codimension d as above. Suppose that x is an 
isolated point of Z n Aa- Then, by [54, Cor. 5.3], there exists a positive integer i 
satisfying that for every sufficiently small Euclidean neighborhood C/ C C" of x 
there is a Euclidean neighborhood F C C'* of a such that for all a' eV 

Aa' is transversal to Z =^ \Z (1 Aa' riU\ = i. (5) 

We call i the intersection multiplicity of Z and Aa at x and wc denote this number 
by i{Z, Aa] x). The multiplicity multa;(Z') of Z at x is defined as the minimum of 
i{Z,Aa;x) over all affine linear subspaces Aa of codimension d such that x is an 
isolated point of Z H Aa [54, Def. 5.9]. It is known that x is a smooth point of Z iff 
mult^(Z) = 1 [54, Cor. 5.15]. 

The following lemma is essential for the first order characterization we are seek- 
ing. 

Lemma 5.8 Let Z C be an algebraic set of dimension d and Aa C C" be 
an aiRne subspace of codimension d, parametrized as above. For x e Z r\ Aa the 
following two conditions are equivalent: 

(a) Aa is transversal to Z at x. 

(b) For every sufRciently small Euclidean neighborhood U C C" of x there is a 

Euclidean neighborhood V C.C^ of a such that for all a' E V the intersection 
Z n Aa' n U contains exactly one point. 

®This is closely related to the question of the expressive power of query languages for constraint 
spatial databases [46]. 

^Mumford considers projective varieties, but the following local considerations clearly hold in 
the affine setting as well. 
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Proof, (a) ^ (b). Assume that ipi{x') = 0, . . . ,ipn-d{x') = are local equations 

of Z at X (i.e., they generate the vanishing ideal of Z in the localization at x). Let 
gi{a,x') = 0, . . . ,gd{a,x') = be equations for Aa, parametrized by the coefficient 
vector a G C^. The transversality of Z and Aa at x implies that the Jacobian matrix 
at x of the polynomial map 

C" C",a;' ^ {(pi{x'), . . .,ipn-d{x'),gi{a,x'),.. .,gd{a,x')) 

is invertible. The implicit function theorem tells us that there is a continuous map 
s : Vb — > C/q between Euclidean open neighborhoods Vq of a and Uq of x such that 
for all a' G Vb, s(a') is the unique solution in Uq of the system of equations 

ipi{x') = 0,...,ipn-d{x') = 0, gi{a,x') = 0, . . . ,gd{a,x') = 0. 

For any Euclidean neighborhood U Q Uq of x, the Euclidean neighborhood V := 
s~^{U) satisfies the statement of condition (b). 

(b) =^ (a). By contraposition, we assume that Aa is not transversal to Z at x 
and show that condition (b) is not satisfied by considering several cases. 

Suppose first that dim^; Z < d. Let U' denote the open neighborhood of x 
consisting of the set of points in C", which do not lie in an irreducible component 
of Z of dimension d. Then Z n A^/ n C/' = for Zariski almost all a' G C'*. If 
condition (b) were satisfied, there would exist sequences Xi ^ x and ai ^ a such 
that e Z n Aa^ and Z n fl ?7' = for all i. Hence Xi U' for all i, which 
contradicts the fact that Xi converges to x. Thus (b) is violated. 

In the following we assume that dim^, Z = d. We may assume that x is an isolated 
point of Z f] Aa since otherwise, (b) is clearly not satisfied. We will distinguishing 
several cases and prove that condition (b) is violated by showing the following claim 
in each case: 

There are two sequences (xi) and (x'-) in C", both converging to x, and 
there is a sequence (oj) in converging to a such that Xj, G Z n Aai (6) 
and Xi 7^ x'^ for all i. 

Let Zi be an irreducible component of Z passing through x such that dim Zi = d. 
If a; is a singular point of Zi, then i{Zi,Aa; x) > nmltx{Zi) > 2 and claim (6) follows 
by the characterization (5) of the multiplicity. 

We may therefore assume that x is a smooth point of Zi. If Aa is not transversal 
to Zi at X, then T^Zi n T^A ^ and therefore T^A contains a line £ tangent to Zi 
at X. There is a sequence of points Xi e Zi, Xi ^ x, converging to x in the Euclidean 
topology such that the secant Sj through x and Xi converges to i. Take A^. to be 
the affine space of codimension d spanned by and Si (here is the orthogonal 
complement of £ in Aa). Since Aq. Z D {x,Xi}, and we can achieve that Oj ^ a 
for a suitable choice of the parameter Oj, the claim (6) follows. 

We are left with the case where Aa is transversal to Zi at x. Since Aa is not 
transversal to Z at x, there must be at least one further irreducible component Z2 
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of Z passing through x. Consider a sequence of points G Z2 — Z\ converging to x 
and such that Xj 7^ x. Consider also points zi, . . . , z^-d in such that the vectors 

linearly independent. Now let Aa^ be the affine space of 
codimension d passing through x^, zi, . . . , z^-d- We can achieve that a. 

On the other hand, since is transversal to Zi at x, we may apply condition (b) 
to Zi and Aa- Passing over to a subsequence of {AaJ, we obtain that there is a 
sequence x'^ G ZiH Aa^ converging to x. This shows the claim (6) and completes the 
proof of the lemma. □ 

In the following, we parametrize a system fi, ■ ■ ■ , fr of polynomials over C by 
its vector of non-zero coefHcients u eC^, and we denote the corresponding zero set 
by Zy_. (Hence we use the sparse encoding, cf. §2.2.) Recall that we parametrize 
affine subspaces Aa C C" of codimension d by elements a E C^. 

Lemma 5.9 For allO < d <n there is a first order formula Fd{u, a) in in prenex 
form with seven quantifier blocks, 0{'n?) bounded variables, and with 0{q + n) 
atomic predicates given by integer polynomials of degree at most 5 and bit-size 
0{1), such that for aU u G C ~ M^g ^j^h dime Zy, = d and aU a G C'': 

Fd{u, a) is true Aa is transversal to Z^- 

Proof. In what follows, we interpret all occur ing formulas over C as first order 
formulas in J^r by encoding a complex number by its real and imaginary part. 

Suppose that A^ is of codimension d. Then property (b) in Lemma 5.8 expressing 
transversality of Aa to Z^ at x can be written as the following formula (^(u, a, x): 

3eo > VO < e < eo 35 > Va' G C'^ 3y G Vz G (||a-a'|| < (5 A 

||y - xll < e A ||z - xll < e A y G n Aa' A {z E Z^^ Aa' =^ y = z)^. 

The property that Aa is transversal to Z„ at all affine intersection points x e ZuCiAa 
then reads as: 

Vx G {x e ZuCi Aa =^ (fiu, a, x)). 
The property that Zu and Aa have no intersection points at infinity is expressed by 

Vx G C"+^ {x eZuAx eAa^xo^O), 

where the bar denotes projective closure (we have now an additional homogenizing 
variable xq). We express the predicate x E Z^ in the form 

Ve>0 3x'GC" 3AgC-{0} (x' G Z„ A ||x - A(1,x')|| < e), 

using the fact that the Zariski-closure of constructible sets equals the Euclidean 
closure. 
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Finally, we can express that codimAa = d by requiring that there exists a linear 
subspace L with dim L > d and A}™ fl L = 0, where denotes the linear space 
associated with Aa- 

Altogether, we see that the transversality condition can be expressed by a for- 
mula in of the required description size. □ 

Remark 5.10 (i) It is not clear whether transversality can be expressed by short 
first order formulas over C since the Euclidean topology is involved. We will 
circumvent this difficulty by working with the first order theory over the reals. 
The next lemma provides a concise first order (over the reals) characterization 
of transversality. However, it is important to keep in mind that we will resort 
to the reals only as a way of reasoning. All computations in the proof of 
Theorem 5.2 will be done by machines over C. 

(ii) Note that the projective closure Z„ is included in but may not be equal to the 
zero set of the homogenization of the polynomials defining Z^- 

5.3 Proof of Theorem 5.2 

We begin with the membership of Degree to FP*^^. Let p = 2q. Then, by 
Theorem 4.9 and Remark 4.10(i), a partial witness sequence a = (q:i, . . . ,a2p+i) 
for the formula Fd{u,a) in Lemma 5.9 can be computed by a machine over C, 
given input {p,k,n,w,m,6,£), in time (nq)'-^^^^ log (5. Note that this quantity is 
polynomially bounded in the sparse input size 0{nqlog5). 

We claim the correctness of the following algorithm for Degree. 

input fi,. . . ,fr with coefficient vector u 

compute d := dim Zu by oracle calls to HNc using Theorem 2.6 
compute a partial witness sequence a = (cki, . . . , oc2p+i) of F^^u, a) 
for i = 1 to 2p + 1 

compute A^j := |Z„ n ^1^^! by an oracle call to #HNc 
compute the majority N of the numbers Ni,. . . , N^p+i 
return N 

Put / := {i G [2p + 1] I Fd{u,OLi) holds}. Lemma 5.9 and Part (ii) of Lemma 5.6 
imply that Aj = deg Z^ for all z G /. Part (i) of Lemma 5.6 tells us that V*a Fci{u, a). 
Since a is a partial witness sequence, this implies that |/| > p (cf. Definition (4.7)). 
This proves the claim. 

It is obvious that the above algorithm can be implemented as a polynomial time 
oracle Turing machine over C. This shows the membership. 

To prove the hardness, note that, by Theorem 3.4, ^^HNc is ^^Pf^-complete. It 
is therefore sufficient to Turing reduce #HNc to Degree. The following reduction 
does so. For a given system of equations first decide whether its solution set Z 
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is zero-dimensional by a call to HNc using Theorem 2.6. This call to HNc can 

be replaced by a call to DEGREE since HNc reduces to DEGREE (recall Z = $ iS 
dcgZ = 0). If dimZ = 0, then compute N := degZ by a call to Degree and 
return N, otherwise return oo. □ 



6 Preliminaries from algebraic and differential topology 

6.1 Euler chctracteristic of compact semialgebraic sets 

It is well known that any compact semialgebraic set S can be triangulated [10, §9.2]. 
Instead of working with triangulations, we will use the more general notion of finite 
cell complexes, since this is necessary for the application of Morse theory in §6.5. 
Compact semialgebraic sets are homeomorphic to finite cell complexes and their 
topology can be studied through the combinatorics of cell complexes. 

We briefly recall the definition of a finite cell complex (also called finite CW- 
complex), see, for instance, [32] for more details. We denote by the closed unit 
ball in M", and by S"'~^ = dD^ its boundary, the (n — l)-dimensional unit sphere. 
An n-disk is a space homemorphic to D"^. By an open n-cell we understand a space 
e" homeomorphic to the open unit ball — dD"^. A (finite) cell complex X is 
obtained by the following inductive procedure. 

We start with a finite discrete set whose points are regarded as 0-cells. 
Inductively, we form the n-skeleton X"' from X^^^ by attaching a finite number of 
open n-cells via continuous maps ipa '■ S""^"*^ This means that X" is the 

quotient space of the disjoint union Uq, D'^ of X""^ with a finite collection of 

n-disks under the identifications x = ipa{x) for x G dD^ = S"-"^. Thus as a set, 
X" = Uq, e", where each is an open n-ccll. We stop this procedure after 

finitely many steps obtaining the compact space X = X'^ of dimension d. 

We note that each cell has a characteristic map : — > X which extends 
the attaching map (pa and is a homeomorphism from the interior of onto e^. 
Namely, we can take to be the composition D^, ^ X""^ U« ^ X"" '-^ X, 
where the middle map is the quotient map defining X". 

Example 1 (i) The n-sphere can be realized as a cell complex with two cells, of 
dimension and n, respectively. The cell e" is attached to e° by the constant 
map ip : S"'~^ e°. 

(ii) Real projective space P"^(M) is defined as the space of all lines through the origin 
in M**"*"^. This is equivalent to identify antipodal points in S"' C M'^"'"^, a pre- 
sentation which in addition yields a natural topology in P'*(M) — the quotient 
topology induced by the identification. Removing the southern hemisphere, 
this is yet equivalent to the space obtained by keeping the northern hemi- 
sphere and identifying antipodal points in the equator. Since the northern 
hemisphere (without the equator) is homeomorphic to e" and the equator 
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with identified antipodal points is just P"~-'^(R), it follows that P"'(M) is ob- 
tained from the n + 1 cells e^, e^, . . . , e" by taking Xq = and, inductively, 
obtaining = P'^(M) from Xk^i by attaching via the identification of 
antipodal points : dD^ X^~^. 

(iii) Complex projective space P"'(C) (already seen in §2.2) is the quotient of the 
unit sphere 5'^"''''^ C C"^^ for the equivalence relation v = \v for all A G C with 
|A| = 1. A reasoning as the one above (taking into account that equivalence 
classes are now homeomorphic to S^) shows that P"'(C) is obtained from the 
n + 1 cells e°, e^, . . . , e^" as above, now getting X2k = P'^(C), for = 0, . . . , n. 

The Euler characteristic of a cell complex X is defined as x{^) = ^'k=o{^^)^^ki 
where is the number of A;-cells of the complex. It is a well-known fact that x{^) 
depends only on the topological space X and not on the cellular decomposition. 
That is, if two cell complexes are homeomorphic, then their Euler characteristics 
are the same. Actually x is even a homotopy invariant. 

Example 1 (continued) For the spaces considered above we obtain, using their 
cell decompositions, that 



A continuous map p: X Y between topological spaces is called a covering 
map if there exists an open cover {Ua} of Y such that for each a, p~^{Ua) is a 
disjoint union of open sets in X, each of which is mapped by p homeomorphically 
onto Ua (see e.g., [12, III. 3]). If the cardinality of the fibre p~^{y) is constant for 
y €Y, then this cardinality is called the number of sheets of the covering map. This 
condition is satisfied when Y is connected. 

Lemma 6.1 If X ^ Y is a covering map with m sheets (m Gnite) and Y is a cell 
complex, then X is also a cell complex and xi^) = i^xO^)- 

Proof. The characteristic maps _D" Y lift to X in exactly m ways. This way, one 
obtains a cell decomposition of X with the number of fc-cells exactly m times the 
number of /c-cells of Y. Thus the alternating sum of these numbers for X is m times 
the alternating sum of these numbers for Y. For more details see [12, Prop. 13.5, 



6.2 Non-compact semialgebraic sets 

There arc several ways to extend the definition of x to non-compact sets. The usual 
one uses singular homology and preserves the property of x of being homotopy 





and x(P"(C)) 



n + 1. 



p. 216]. 



□ 
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invariant. In §6.3 we will see another way which does not, but instead has a useful 
additivity property. 

In algebraic topology one assigns to a topological space X and a field F the 
singular homology vector spaces Hk{X\F) for A: € N, which depend only on the 
homotopy type of X and F. The kth Betti number over F bk{X; F) of X is defined 
as the dimension of Hk{X;F). In case F = Q wc write bk{X) and talk about the 
kth Betti number of X. The Euler characteristic of the space X is defined by 

xiX) = J2i-'^)''dimpHkiX;F) (7) 

fceN 

(if this sum is finite). The Betti numbers bk{X]F) depend on the field F as well 
as on X. Remarkably, their alternate sum, is independent of F. In addition, for 
cell complexes X, this alternate sum coincides with x{^) as defined in §6.1. For a 
general reference to homology we refer to [32, 55]. 

More generally, one can assign to a pair Y <Z X oi topological spaces the 
relative Euler characteristic xi^-,'^) •= x(^) ~ x(^)- It can also be character- 
ized in terms of the relative homology vector spaces H]^{X,Y; F) as xi^^^) = 
Y^j^^^{—l)''dimpHk{X,Y;F). Since Hk{X,Y; F) depends only on the homotopy 
type of the pair (X, Y), the same holds for the relative Euler characteristic x(X, Y). 
Note that Hk{X,$;F) = Hk{X;F) and x{X,0) = x{X). 

Lemma 6.2 Let Z he a compact real algebraic n-dimensional manifold and K C. Z 
be a compact semialgebraic subset. Then 

^ly _K\ = l ^^^^ ~ ^^^^ 

I X{K) ifn is odd. 

Proof. A fundamental duality principle going back to Poincare and extended by 
Alexander and Lefschetz states that for an n-dimensional manifold Z and a com- 
pact subset K carrying the structure of a cell complex, the relative homology space 
Hf:(Z,Z — K;7j2) is isomorphic to the homology space^ i?„_fc(-fC; Z2) for all k. 
See [32, Prop. 3.46, p. 256] or [12, Thm. 8.3, p. 351]. Therefore, we get under 
the assumptions of the theorem that 

x{Z)-x{Z-K) = x{Z,Z-K) = ^{-l)''dimHk{Z,Z-K;Z2) 

k 

= (-ir5;(-l)^dimif,(A-;Z2) = (-irx(i^)- 
k 

This implies the claim in the case where n is even. When n is odd, we obtain that 

x{Z — K) = xi^) + xiZ). On the other hand, applying the above formula for 
K = Z yields x(Z) = -xiZ) and thus x{Z) = 0. Hence x{Z - K) = x{K). □ 

Actually, one gets a natural isomorphism with the cohomology vector space i7"~'°(X;Z2) in- 
duced by the Z2-orientation of the manifold Z, but this is not important for our purposes. 
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6.3 Modified Euler characteristic 



Let S be the disjoint union of two semialgebraic sets Si and 82- In general, it is 
not true that x(S') = xi'^i) + x('S'2) • For a counterexample, consider the closed 
3-dimensional unit ball decomposed into its interior and its boundary S'^. 

Yao [73] defined the modified Euler characteristic x* of semi-algebraic sets, which 
satisfies an additivity property, and coincides with the usual Euler characteristic for 
compact semialgebraic sets. The following proposition from [73] characterizes this 
notion. 

Proposition 6.3 There is a unique function x* mapping semialgebraic sets to in- 
tegers, which satisGes the following properties: 

(i) IfS = Uili '^i a disjoint union of semialgebraic sets, then x*(S') = Yli^i X*{'S)- 

(ii) We have x*(<S') = x(-S') for compact semialgebraic sets. 

(iii) If there is a semialgebraic homeomorphism S ^ T, then x*{S) = x*{T)- 

Proof. For the proof of existence, which relies on Hironaka's triangulation theo- 
rem [37] for bounded (not necessarily closed) semialgebraic sets, we refer to [73]. 

The proof of uniqueness shows that, in principle, the computation of x* can be 
reduced to computations of x for compact semialgebraic sets. Since this is useful 
for calculating some examples, and to familiarize the reader with the notion of the 
modified Euler characteristic, we present the simple proof of uniqueness. 

Any unbounded semialgebraic set S C R" is semialgebraically homeomorphic to 
a bounded one. Namely, S is homeomorphic to its its image under the inverse of the 
stereographic projection 5" — {(0, ... ,0, 1)} —>■ M",x 1— > y given by the equations 
Vi = Xi/{1 — Xn+i). Therefore, by property (hi), it suffices to show uniqueness for 
bounded semialgebraic sets S. We proceed by induction on the dimension of S. The 
case dimS* < is clear. Consider the disjoint union S = S U R, where R := S — S . 
We have dimR < dimt^S" < dimS" since R is contained in the boundary dS of S, 
cf. [10, Prop. 2.8.12]. Since S is compact we have x*{S) = xi^) by property (h). 
Property (i) implies that X*(<S') = x(S') — x*(-R), hence x*('S') is determined by 
X*(i?), which in turn is uniquely determined by the induction hypothesis. □ 

Example 2 The inverse image of R" under the stereographic projection is S'^ minus 
a point, hence x*(l^"') = x{S^) - 1 = (— l)"- Note that, in contrast with x, X* is 
not invariant under homotopies. 

Corollary 6.4 If Si, . . . , Sn are semialgebraic subsets of M", then we have 



N 




where the summation is over all nonempty subsets I of [N]. 
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Proof. This follows from the inclusion-exclusion principle taking into account that 
X* behaves additively with respect to disjoint unions. □ 



6.4 Locally closed spaces and Borel-Moore homology 

A noncompact locally closed set S can be compactified by adding just one point. 
More specifically, there is a compact semi-algcbraic set S and a continuous semi- 
algebraic map l: S ^ S, which is a homeomorphism onto its image, such that 
S — i{S) consists of just one point oo, cf. [10, 2.5.9]. 

Let 5 be a locally closed semialgebraic set and F be a field. If S is not compact, 
then the Borel-Moore homology vector spaces of S over F are defined as the relative 
homology spaces of the pair {S,oo), that is, Hf^{S;F) := Hk{S, oo; F), cf. [10, 
§11.4]. If S is compact, then we define i?f^(S';F) = HkiS,F). Moreover, we 
define the fcth Borel-Moore Betti number of S, denoted 6™ {S) , as the dimension of 
-fffe'^(5';Q). Thus we have bl^{S) = bk{S) for compact S. 

Prom the above, the following well-known characterization easily follows. 

Proposition 6.5 Let S be a locally closed semialgebraic set. Then 

x*{s) = ^{-i)'br{s). 

fceN 

Proof. If S is compact the result is trivial. Otherwise, we have x*{S) = X*{S) — 
X*(oo) = x{S) — x{oo) = x{S, oo). On the other hand 

X{S, oo) = ^(-1)^ dimHkiS, oo; Q) = J^(-l)'= dimHr{S; Q), 

fceN fceN 

which shows the assertion. □ 



Remark 6.6 (i) x*{S) can also be interpreted as the Euler characteristic of S with 
respect to the cohomology H*{S;Q) of S with compact supports, a notion 
naturally occuring in the Poincare duality theorem for noncompact manifolds, 
cf. [32, §3.3, p. 242). 

(ii) It is an important fact that for a complex algebraic variety W we have x*(W^) = 
x(H^). If W is smooth of complex dimension n, then this follows from the 
Poincare duahty Hk{W) ~ H^"'~'^{W), using the interpretation of x*0^) ^ 
the Euler characteristic of the cohomology H*{W) with compact support. For 
the proof of the general case see [25, Exercise §4.5, p. 95 and Notes §4.13, 
p. 141]. 
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6.5 Morse Theory 

We recall now some notions and facts from Morse theory. A general reference for 
this is [51]. 

Let Z he a differentiable manifold and (p: Z ^ W he differentiable. A point 
X & Z is a critical point of (p if the differential dxP- T^Z M vanishes. In this 
case, one may consider the Hessian Hxip: T^Z x T^Z M. of (p at x, which is 
a symmetric bilinear form (defined by the second order derivatives of (p in local 
coordinates). The function (p is called nondegenerate at the critical point x if its 
Hessian is nondegenerate at x. The function (p is called a Morse function if all its 
critical points are nondegenerate. 

We call the number of negative eigenvalues of a symmetric matrix or of a sym- 
metric bilinear form its index. The index of 99 at x is defined as the index of HxP- 
Throughout the paper, we will use the convenient notation {ip < r} := {x & Z \ 
(p{x) < r}. 

The main theorem of Morse theory [51, Thm. 3.5] states the following. 

Theorem 6.7 Assume that ip: Z — >■ M is a Morse function on a diffcrcntiahle mani- 
fold Z with finitely many critical points. Moreover, assume that {ip < r} is compact 
for all r S M. Then Z has the homotopy type of a cell complex with one cell of 
dimension k for each critical point of ip of index k. □ 

We will use the following consequence of this result, adapted to the semialgebraic 
setting. 

Corollary 6.8 Let Z C R" he a real algebraic manifold. Then, 

(i) The Euclidean distance function La'. Z ^ M., x \\x — a|p, is a Morse function 

for Zariski almost all a G M". 

(ii) Suppose that La is a Morse function on Z. Then the number of critical 

points of La with index k is Enite for allQ < k < n and Ylk=oi~^)'' equals 
the Euler characteristic x{^) of Z. 

Proof, (i) The first claim follows as in [51, §6] by using the semialgebraic Morse- 
Sard Theorem [10, Thm. 9.5.2]. 

(ii) It is easy to see that the set of critical points of La is semialgebraic. 
Moreover, critical points are isolated. Since semialgebraic sets have finitely many 
components, it follows that there are only finitely many critical points. Note that 
Zr]{x G M" I La{x) < r} is compact for all r E M. Hence we can apply Theorem 6.7 
and the claim follows from the definition of □ 

Let Jif be the set of polynomials / G M[Xi, . . . ,X„] satisfying that Z{f) / 
along with the regularity condition 

VxGM" (/(a;) = 0^grad/(x)7^0). (8) 
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Note that Z{f) is a smooth hypersurface for / G Jf . 

Consider / G and Z = 2(f)- Then x € Z is a critical point of La if and 
only if dxf^f{x){xk — Ofc) = 0. Let a; be a critical point of La such that (w.l.o.g.) 
dx„f{x) 7^ 0. By the implicit function theorem, locally around x, Z is the graph of 
a function (ti, . . . ,tn-i) ^ uiti, ■ ■ ■ ,tn-i) which defines a local coordinate system 
around x, 

{ti,...,tn-i) ^x{t) := 



Lemma 6.9 The Hessian {Hij) = {^t^^t ■La{x{t))) of the distance function La at x 
in terms of the local coordinates ti is given hy 

1 

{dxjfS,, + {dxj){dx,f) + (xn - an){{dxj){dx,dxj) - {dx.^dxj){dxj)). 
Proof. By differentiating La{x) = '^k{xk — a-kf' with respect to t^ we obtain 

n 

dtiLa = 2 ^(xfe - ak) dtiXk = 2{ti -ai + {y- an)dtiy). 
k=l 

Differentiating again with respect to tj yields 

Hij = dtAj^a = 2((^ij + {dtiy){dt^y) + {y - an) dtAjV)- 

Prom f{ti, . . . ,tn-i,y{ti, . . . ,tn-i)) = we get dt^y = by differentiating. 

Differentiating this again with respect to tj we obtain 

. . _ -{dt,dxj){dxj) + {dxj){dt,dxj) 

By plugging these expressions for the partial derivatives of y into the above formula 
for Hij and taking into account that ti = Xi for i < nwe obtain the asserted formula. 

□ 

As in Section 5, we denote by G the vector of non-zero coefficients of the 
polyomial / = /« of degree S in Xi, . . . and write Z„ := Z(fu) for its zero set 
in M". 

The following lemma gives a certificate for La to be a Morse function on Zu in 
form of a parametrized first order formula. It plays a similar role for the completeness 
proof of EULER^ as the certificate for transversality for the completeness proof of 
Degree, which was provided in Lemma 5.9. 

Lemma 6.10 There is a first order formula F{u, a) in in prenex form with one 

quantifier block, n bounded variables, and with 0{n) atomic predicates given by 
integer polynomials of degree at most 0{n6) and bit-size 0{nlog{np)) such that, 
for allueW such that fu&J^ and all a G M", the following holds: 

F{u, a) is true <^=^ La'. ^ M is a Morse function. 
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Proof. The fact that L^: — >^ M is a Morse function can be expressed by the 
fohowing formula 

(n n \ 

k=l fc=l ^ 

where, we recall, HxLa denotes the Hessian of La at x. We now replace HxLa by the 
explicit expression for it given in Lemma 6.9, after making the appropriate changes 
due to the fact that we require the fcth partial derivative of / to be nonvanishing 
at X instead of the nth derivative. The assertion follows now easily by inspecting 
the above formula. □ 

7 Complexity of the Euler characteristic 

Another main result of this paper proves the completeness in FPj^ * of the following 
problem over M. 

EuLER^ {Modified Euler characteristic) Given a semialgebraic set 5 C R" as a 
union of basic semialgebraic sets 

S = [JUAx e M" I gi{x) = 0, fn{x) > 0, . . . , fir,{x) > 0}, 

decide whether S is empty and if not, compute x*{S)- 

Theorem 7.1 The problem Euler^ is FPj^ ^-complete with respect to Turing re- 
ductions. 

The upper bound in Theorem 7.1 is proved in several steps: in Section 7.1 we 
reduce the basic semialgebraic case to the case of a smooth hypersurface. This 
case is then treated in Section 7.2 based on Morse theory and the concept of partial 
witness sequence developed in Section 4.3. Finally, we combine these two ingredients 
in Section 7.3 to treat the case of arbitrary semialgebraic sets, using the inclusion- 
exclusion principle, which is possible due to the additivity property of the modified 
Euler characteristic. 



7.1 Basic semialgebraic, projective and affine varieties 

Lemma 7.2 Let g, fi, . . . , fr G . . . , Xn] be of degree at most 5 and S := {x & 

M" I g{x) = 0,fi{x) > 0, . . . , fr{x) > 0}. Put go := g and dehne for I < i < r 

r 

9i '■= ^n+ifi ~ 1) '= Xo^^9i{^l/^0, ■ ■ ■ , Xn+r/Xo), H :=^2 ^j- 

i=0 

Then, $ := Z{H — 1) C M""'"''"'"^ is a smooth afhne hypersurface and 

\n+r 
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Proof. Note that, for i = l,...,r, Gi G R[Xo, . . . , Xn^r] is homogeneous and 
gi G . . . ,Xn^r]- Define the affine variety Ya and the projective variety Yp by 

Ya := Z{go, ...,gr)^ Yp := Z{Go, ...,Gr) = Z{H) C P«+'-(R). 

For e G {—1, l}*" consider the open subsets Y^ := Ya n (n[^]^{sgn(a;„+j) = gj}) of y^. 
Clearly, each 1^ is semialgebraically homeomorphic to 5*. Moreover, Ya is the disjoint 
union of the Y^. Hence 

2'x*{S) = Y,X%Y,)=x*{Ya). (9) 

e 

Consider the open subset V := Ypf] {Xq / 0} of Yp, which is semialgebraically 
homeomorphic to Ya. Since we homogenized with exponent (5 + 3, which is one 
higher than the maximum degree 5 + 2 of the gi, wc have Yp — V = 2^pn+r(]j)(Xo) — 
P"+''-i(M). By additivity of x* we have x{Yp) = x(I"*+''"^ W) + X*{V), hence 

y*^v ) = = f if n + r is even, 

^ X It^J ifn + risodd. ^ ^ 

Note that 1 is a regular value of H, since H = {deg H)Y^^XidxiH by the 
homogeneity of H. Hence the "Milnor fibre" 

$ := {x G R^+'^+i I H{x) = 1} 

is a smooth affine hypersurface. Put U := {x e P"'+''(M) | H{x) ^ 0}. We claim 
that the canonical map 

TT: ^ ^U, [xq,..., Xn+r) ^ {xq : ■ ■ ■ : Xn+r) 

is a covering map with two sheets. Indeed, Tr~^{U fl {Xi ^ 0}) = ($ fl {Xi > 
0}) U ($ n {Xi < 0}), and tt induces homcomorphisms from both $ fl {Xi > 0} and 
$ n {Xj < 0} to [/ n {Xj / 0}, respectively. 

By Lemma 6.1 we have x(^) = '^xiU). On the other hand, by Lemma 6.2 and 
Example 1, we get x{U) = 1 — x{Yp) if n + r is even and x(C/) = x(^) if n + r is 
odd. Altogether, we obtain 

^(Y \ - S - ifn + riseven, 

1 5X($) ifn + risodd. ^ ' 

Combining Equations (9), (10), and (11) the assertion follows. □ 
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7.2 The c£ise of a smooth real hypersurface 

Consider the function xjt : ^ — / i-^ x{^{f)) computing the Euler character- 
istic of the smooth hypersurface Z{f) given by / G Jf. Note that we don't consider 
the modified Euler characteristic here. 

Proposition 7.3 The function xjr belongs to FP^ 

Proof. Let INDEX be the following decision problem. An input to INDEX is a 
tuple (u,a,x,k, J), where u encodes a real polynomial f inn variables, a,x € M", 
k e N and J C [n] is nonempty. The question is to decide whether x is a critical 
point of index k of the function La : Z.^ ^ M satisfying dxjfix) 7^ for all j G J. 

The problem INDEX is in Pjk. Indeed, given the tuple {u,a,x,k, J), one first 
computes the Hessian H^La by using Lemma 6.9. Then, one computes its char- 
acteristic polynomial (a computation known to be in FPr, see [8, 16], and finally 
one uses Sturm's algorithm to compute the number of real zeros in the interval 
(—00,0) (again in FPr, see [26]). Comparing this number with k decides INDEX 
for (u, a, X, k, J). 

Given {u, a), let x+{u, a) denote the number of {x, k, J) such that {u, a, x, k, J) G 
INDEX and k + |J| is odd. Similarly, we define X-(''^)0) by requiring that k + |J| 
is even. Since INDEX G Pr, the functions M°° x R°° ^ N U {00} mapping {u, a) to 
X+{u,a) and X-(,u,a), respectively, are in ^P^. 

We claim that 

xi^u) = X+i'^j 0^) ~ X-i^j ^) if is a Morse function on Z^- (12) 

In fact, if Nk denotes the number of critical points of La on Z^, we have xi^u) = 
^^( — l)^A^fe by Corollary 6.8. Let N^j denote the number of critical points x of La 
on Zu satisfying dxjf{x) 7^ for all j G J . By the inclusion-exclusion principle we 
have Nk = Z^j^0(— l)''^'~^-^fcj (note that at every point of Zu at least one partial 
derivative does not vanish due to the regularity condition for /„). This implies 

X(Z„) = ^(-l)'=iVfc = (-l)'+l'"-'iVfcJ = X+(^,a)-X-(n,a) 

k k,J^9 

as claimed in (12). 

Lemma 6.10 and Theorem 4.9 imply that a partial witness sequence a for the 
first order formula F{u,a) certifying that L^: ^ R is a Morse function can be 
compTitcd (uniformly) by a division-free straight-line program with (np)'^(i) log((5) 
arithmetic operations, using 1 as the only constant. 

The following algorithm computing XM' can be implemented as a polynomial 
time oracle Turing machine querying oracles in #P]r. 
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input / G encoded by its coefficient vector u 

compute a partial witness sequence o: = {cxi, . . . , 0:2^+1) of F{u, a) 

for £ = 1 to 2p + 1 

compute x{u, oce) := x+{u, oie) - X-{u, Oie) 
compute the majority x('u) of the numbers x(^) cki); ■ ■ ■ ) x('f^) Q!2p+i) 
return x(u) 

In order to show that this algorithm actually computes the Euler characteristic 
of its input, put A := {i (z [2p+l] | F{u, cxf ) holds}. By definition of F we know that 
Lot( is a Morse function on Z„ for all £ € A. Hence, by (12), xi^u) = ct^) for all 
£ e A. On the other hand, by Proposition 6.8(i) we have V*aF(ii, a). Since a is a 
partial witness sequence, this implies that |A| > p (cf. Definition (4.7)). Therefore, 
the algorithm indeed computes the Euler characteristic oi Z^- □ 

7.3 Arbitrary semialgebraic sets 

Proposition 7.4 The problem EULER^ is contained in FP^^*. 

Proof. Consider an instance S = U^^^Si of the problem Euler^, where t > 1 and 
Si = {x e \ gi{x) = 0, fii{x) > 0, . . . , finix) > 0}. Emptyness of S can be easily 
decided in FP*^«. 

Assume S" 7^ 0. By adding dummy inequalities 1 > 0, we may assume that = r 
for all i. Corollary 6.4 tells us that 



We will assume that each Sj is described by exactly rt inequalities, which can be 
achieved by adding further dummy inequalities 

According to Lemma 7.2, we can assign to each nonempty index set / C [t] a 
homogeneous polynomial Hj G M[Xo, . . . , Xn+rt], such that x*{Si) can be expressed 
by the Euler characteristic of the smooth affine hypersurface $7 := Z{Hj — 1) in 
j^„+l+rt follows 



x*(5) = E7^0(-i)i'i-V(5/) 



(13) 



where for / C [t], the basic semialgebraic set Sj C is defined by 





(14) 




(15) 
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We proceed now similarly as in the proof of Proposition 7.3. Let p be the number 

of real parameters of all the polynomials gi, fij involved in the above description of 
the set S. To emphasize the dependence on u, we will write <^/^„ instead of For 
a projection point a G ^^d a parameter u € W we consider the distance 

function La : — > M, x h- > ||x — ap. 

Similarly as in the proof of Proposition 7.3, we assign to n G W, a G R""^^"*'''*, 
and / C [t] two values x+,iiu,a),X-,i{u,a) G N such that (cf. (12)) 

X{^i,u) = X+,i{u, a) - x-,i{u, a) if La is a Morse function on (16) 

Namely, x+,l{u, a) is defined as the number of triples {x, k, J), where x G R"+^+''*, 
A; G N, 7^ J C [n + 1 + rt] such that x is a critical point of index k of the function 
La - ^i,u R satisfying dxjf{x) 7^ for all j (z J and such that k + |J| is odd. 
Similarly, one defines X-,i(,u,a) by requiring that k + \ J\ is even. As in the proof 
of Proposition 7.3, one shows that the functions {0, 1}°° x M°° x M°° ^ N U {00} 
mapping {I,u,a) to x+,i{u,ci) and X-,i{u,a), respectively, are in #Pr. 

Assume now that n, a are chosen such that La is a Morse function on ^i^u for 
all nonempty subsets / of [t]. Plugging (16) into (15) we obtain 

(-ir+^*2^*+ix*(^) = 2+j;(-l)l^l(x+,/(ix,a)-X-,/(^x,a)) = 2+x+{u,a)-xA^,a), 

where we have put 

X+{u,a):= x+,i{u,a), %-('",«):= Yl X-,i{u,a). 

I^ID,\I\ even //0,|/| odd 

According to Lemma 3.9, the functions {u,a) h- > x+(w, a) and {u,a) 1— > x-{u,a) are 
in #Pm. 

Consider the first order formula Gi{u, a) in provided by Lemma 6.10, which 
expresses the fact that La'. ^ R is a Morse function. Define the first order 
formula G{u,a) := Aj^ijfGi{u, a), which certifies that, for all nonempty index sets 
/ C [t], La'. — > M is a Morse function. Theorem 4.9 and Remark 4.10 imply 
that a partial witness sequence a = (cti, . . . , Q!2p+i) for the formula G{u, a) can be 
computed in time polynomial in the input size of S. (Note that it does not harm 
that the number of atomic predicates of G{u, a) is exponential in the input size 
of S.) 

After all these preparations, we see that the modified Euler characteristic of S 
can be computed by essentially the same algorithm as in the proof of Proposi- 
tion 7.3. The modifications are as follows: replace the formula F by G, reinterpret 
X+(w, a), X+('Ji! o) in the above way, and return (— 1)""'"''*2~^*~-^(2 + x(it)) where, 
again, x('") is obtained by taJiing a majority vote on the x+{u,Oii) — X-{u,cXi). 
This algorithm can be implemented as a polynomial time oracle Turing machine 
accessing oracles in #P]k. The proof of correctness is identical as for the proof of 
Proposition 7.3. □ 
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Proof of Theorem 7.1. The membership of EULER^ to FPJ * is the content of 
Proposition 7.4. By Theorem 3.4, t^FeaSr is ^^Pj^-complctc. To prove the Turing- 
hardness of EuLERj^ for T^P]K, it is therefore sufficient to Turing reduce ^Feas^ to 
EulerJ. The following reduction does so. For a given real polynomial first decide 
whether its solution set Z is zero-dimensional by a call to FeaSr using Theorem 2.6. 
This call to FeaSr can be replaced by a call to EuLERj^ since FeaSk reduces to 
EuLERj^ (this follows from the case distinction in the definition of the problem 
EuLER^). If dimZ = 0, then compute N := x*{Z) by a cah to EULER^ and 
return N, otherwise return oo. □ 

Remark 7.5 In the papers [13, 66], the Euler characteristic of a real algebraic va- 
riety is expressed by the index of an associated gradient vector field at zero, which 
can be algebraically computed according to [23]. Although Morse theory is not ex- 
plicitly mentioned in [13, 66], the main idea behind these papers is an application 
of this theory as exposed in [53]. The single exponential time algorithm in [3] for 
computing the Euler characteristic uses Morse theory explicitly and in a crucial way. 
However, we note that the reduction in [3] from the case of an arbitrary semialge- 
braic set to the case of a smooth hypersurface, as well as the reductions in [13, 66], 
cannot be used in our context, since it is not clear how to compute the deforma- 
tion parameter or the sufficiently small radius of the intersecting sphere within the 
allowed resources (polynomial time for real machines). Instead, we have expressed 
the Euler characteristic of a real projective variety by the Euler characteristic of its 
complement, which in turn can be expressed as the Euler characteristic of a "Milnor 
fibre" , which is a smooth hypersurface. 

8 Completeness results in the Turing model 

It is common to restrict the input polynomials in the problems considered so far to 
polynomials with integer coefficients. The resulting problems can be encoded in a 
finite alphabet and studied in the classical Turing setting. In general, if L denotes a 
problem defined over R or C, we denote its restriction to integer inputs by L'^. This 
way, the discrete problems HN^, DlM^, Degree°, Euler^°, etc. are well defined. 

We are going to show next that all the above problems are (Turing-) complete in 
certain discrete complexity classes. These classes are obtained from real or complex 
complexity classes by the operation of taking the Boolean part. 

8.1 Basic complete problems in Boolean parts 

A problem that has attracted much attention in real (or complex) complexity is 
the computation of Boolean parts [15, 20, 21, 22, 39, 43]. Roughly speaking, this 
amounts to characterize, in terms of classical complexity classes, the power of re- 
source bounded machines over M or C when their inputs are restricted to be binary. 
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Definition 8.1 Let C be a complexity class of decision problems over M or C. Its 
Boolean part is the classical complexity class 

BP(C) := {5n{0,l}°° \S eC}. 

The study of Boolean parts has been successful in the setting of additive ma- 
chines, where practically all natural complexity classes have had their Boolean parts 
characterized [17, 22, 39]. In contrast, much less is known in the setting of unre- 
stricted machines. Two of the most significant results state that BP(Pc) C P'^'^ [21] 
and BP(PAR]r) = PSPACE/poiy [20], and a third one is discussed in Proposition 8.3 
below. For stating it, recall that RP denotes the classical complexity class of prob- 
lems decidable by randomized machines in polynomial time with (one-sided) error. 
It is well-known that P'^^ C n^, where 11^ denotes a class in the second level of the 
polynomial hierarchy (see [2, 58] for details). 

The following upper bound for HNc° was obtained by Koiran [40] . 

Theorem 8.2 HNc'^ belongs to RP^^ (and therefore to iP) under the generahzed 
Riemann hypothesis GRH. □ 

A natural restriction for real or complex machines (considered e.g. in [22, 39, 43]) 
is the requirement that no constants other than and 1 appear in the machine 
program. Complexity classes arising by considering such constant-free machines are 
indicated by a superscript as in F^, NP]^, etc. Note that wc have already used a 
superscript for a decision or counting problem L over M or C, denoting by its 
corresponding problem restricted to integer inputs. Thus, the superscript means 
restriction to bits in two different contexts: machine constants when considering 
complexity classes, and input data when considering problems. This should not 
create any confusion. 

Theorem 8.2 provides an upper bound for HNc°. On the other hand, the clear 
NP-hardncss of HNc'^ provides a lower bound. Yet there is a gap between NP and 
RP^*^ and the problem of how to close it (with regard to HNc*') is, as of today, an 
open question. The following result elaborates on that question. 

Proposition 8.3 (i) HNc° and DlMc° are BP(NP^)-comp]ete. 
(ii) Assuming GRH, we have NP C BP(NP^) C RP^"". 

Proof. The completeness of HNc° in part (i) follows from the following fact. The 
FPc-reduction from an arbitrary NPc-problem to HNc exhibited in [8], when ap- 
plied to a problem L in NP^?, yields a FP-reduction from to HNc*^. This shows 
that HNc'^ is BP(NP£)-complete. The completeness of DlMc'^ follows from Theo- 
rem 2.6(i). 

For the reasoning above to hold it is essential that we only consider problems 

defined by NPc-machines that do not use complex constants. Otherwise, these 
constants would appear as coefficients in the constructed polynomial system. 
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The second inclusion in part (ii) follows from part (i) and Theorem 8.2. The 
first inclusion is trivial. □ 

It is believed [58, p. 255] that RP'^^ has no complete problems. Thus, it follows 
from Proposition 8.3 that the equality BP(NP^) = RP^"^ is unlikely to hold. 

The rest of this section is devoted to completeness results in Boolean parts in 
the spirit of Proposition 8.3. Before stating our result, we note that the definition 
of the Boolean part can be extended to classes such as #Pc or ^P^ in an obvious 
way. Thus, we define the class of geometric counting complex problems as GCC := 
BP(#P^) and the class of geometric counting real problems GCR := BP(#P]g'). These 
are classes of discrete counting problems, closed under parsimonius reductions, which 
can be located in a small region in the general landscape of classical complexity 
classes. Namely, we have 

#P C GCC C GCR C FPSPACE, 

where the rightmost inclusion follows from Theorem 3.7 and [20]. 

Proposition 8.4 (i) Feas^, SASm°, and DlM^ are BP(NP^)-I^rmg-compJete. 

(ii) #SASr° and #Fea^ are GC R-T^rmg-compJete. 

(iii) ^HNc" is GCC-Turing-complete. 

Proof. For the hardness in part (i) we use the argument in the proof of Proposi- 
tion 8.3(i), namely, that the reductions from an arbitrary NPjj-problem to FeaSr 
or SASr yield reductions from problems in BP(NP^) to FeaSr^ or SASr'^, respec- 
tively. For the hardness in parts (ii) and (iii) one uses the reductions in the proof 
of Theorem 3.4. The memberships in all statements are clear except for DlM^, for 
which the claim follows from Theorem 2.6(ii). □ 

Remark 8.5 One can show that BP(NP^) = BP(NPc) and GCC = BP(#P^) = 
BP(#P£;). Hence it is immaterial whether we allow the use of complex machine 

constants in the definition of these classes or not. Moreover, it is possible to extend 
Proposition 8.3(fi) to BP(FP^^'^) C RP^^, assuming GRH. The proof relics on the 
possibility to eliminate complex constants using witness sequences, as developed in 
[7, 41, 44]. Details will be given elsewhere. 

We can give some evidence that counting over C is indeed harder than deciding 
feasibility over C. 

Corollary 8.6 If #Pc C FP^^^, then the classical polynomial hierarchy collapses 
at the second level, assuming GRH. 
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Proof. Assuming C FP^^ and taking Boolean parts, we get by Remark 8.5 
#P C BP(#Pc) C BP(FPc*) C Rpl^P C tf. 

Toda's theorem [68] states that PH C P#'^. Hence we conclude PH = 11^, which 
means that the polynomial hierachy collapses at the second level. □ 

8.2 Degree and Euler characteristic in the Turing model 

We can now easily deduce completeness results for the discrete versions of the prob- 
lems to compute the degree or the modified Euler characteristic. 

Theorem 8.7 (i) Degree° is FP^'^^-compiete with respect to Taring reductions. 
(ii) EULER^° is FP^^^ -complete with respect to Turing reductions. 

Proof, (i) The proof given in Section 5 for the membership of DEGREE to FP^^"- 
applies in our case with only one modification. The algorithm in the proof of The- 
orem 5.2 computes the partial witness sequence ot (this is done in FPc) and then 
performs 2p+l oracle calls to #Pc to obtain the numbers Ni for i G [2p + 1] . While 
it is clear that the computation of a is in BP (FPc), it is equally clear that it is not 
in FP due to the exponential coefficient growth caused by repeated powering (cf. 
Lemma 4.5). A way to solve this is to "move" the computation of a to the query. 
That is, one considers the problem of computing iVj with input {u,i). Clearly, this 
problem is in BP(7^Pc): one first computes a in FPc and then Ni in t^P^- 

The hardness of Degree*^ follows as in Theorem 5.2 using the second statement 
in Theorem 2.6(i) instead of the first. 

(ii) The proof for EULER^° is a modification of the proof of Theorem 7.1, similar 
as for part (i). □ 



Remark 8.8 (i) The algorithms for Degree and EuLER^ above can be further 
simplified. Since we can bound the description size of the formula F{u,a) 
or G{u, a) by taking into account a bound on the bit-size of the components 
of the given u G Z^', the input vector u does not need to be considered as 
a parameter any more. Therefore, we may take p = 0. The partial witness 
sequence then consists of a single vector a G Z'^ and only one oracle call to 
#HNc° (or two oracle calls to #Fea^) are needed. 

(ii) Alternatively, the algorithm in the proof of Theorem 5.2 (or Theorem 7.1) could 
be modified as follows. By part (i) we may assume that p = 0. The straight- 
line computation for the partial witness a £ of F cannot be executed in the 
bit model because of the exponential coefficient growth. However, we can easily 
remedy this by describing the construction of the partial witness sequence by 
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existentially quantifying over additional variables f3i, . . . ,f3q along the recursive 
description in Lemma 4.5. We then query #HNc° for the system of equations 
in the variables x, on and (3i, . . . ,Pq expressing the recursive construction of 
Qfi and the fact that x G D L^. 

In the Turing model we can also prove a completeness result for the computation 
of the (non-modified) Euler characteristic: consider the problem 

EULER]g (Euler characteristic for basic semialgehraic sets) Given a basic semial- 
gebraic set 5 = {re G M"" | g{x) = 0, fi{x) > 0, . . . , fr{x) > 0}, decide whether 
S is empty and if not, compute xi'^)- 

Theorem 8.9 Euler^ is FP^^'^-compiete with respect to Turing reductions. 

To prepare for the proof, recall that a closed scmialgcbraic set S C M" has a 
conic structure at infinity [5, Prop. 5.50], which implies that there exists r > 
such that for all r' > r there is a semialgehraic deformation retraction from S to 
Sr' := Sr\{x € M" I ||x|| < r'}. We will call r a cone radius of S at infinity. Clearly, 
we have x{S) = x{Sr) = X*{Sr)- 

Lemma 8.10 Let p G Z[Xi, . . . ,Xn] be of degree at most 6 with coefficients of 
hit-size at most £. Then, there exist m = (n5t)^^^^ such that 2^"" is a cone radius 
of Z{p) at infinity. 

Proof. (Sketch) In [27] it is shown that there is a first order formula $(r) in in 
prenex form with the free variable r such that there exists ro > with 

[ro, oo[C {r G M I $(r) true} C {r € M | r is a cone radius of Z{f) at infinity}. 

By an inspection of the constructions in [27, 59] one can show that the formula $(r) 
has a bounded number of quantifier blocks, n^^^^ bounded variables, and m atomic 
predicates given by integer polynomials of degree at most d and bit-size at most 
I' such that \og{dml') < {n5i)^^^\ The tedious details of verifying this statement 
about the desription size of $(r) are omitted for lack of space and left to the reader. 

According to Theorem 4.1, the formula -i$(r) is equivalent to a quantifier-free 
formula in disjunctive normal form v|^^ ^j=i {hij{r)Aij(}), containing integer poly- 
nomials hij{r) of bit-size at most L such that logL < {n5i)'-^^^\ 

Let p € M be the maximum of the real roots of the nonzero hij. We have 
P < 1 + ll^lloo < 1 + 2^. Note that the sign of hij{x) is constant for x > p. Therefore, 
since the set {r > | -i<I>(r)} is bounded, we have {r > | -'$(r)} C]0, p]. Hence 
2 -I- 2^ is a cone radius of Z{f) at infinity, which proves the claim. □ 

Proof of Theorem 8.9. The hardness of EuLER^ follows as in the proof of Theo- 
rem 7.1. We prove now that Euler^ belongs to FP^*''^. For given 5 = {x € M" | 
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g{x) = 0, > 0, . . . , fr{x) > 0}, compute the polynomial 
p{X,Y) ■.= g{Xf + j2{y'fi{X)-lf 

i=l 

in the variables Xi, . . . , Xn, Yi, . . . , Y^. As in the proof of Lemma 7.2 we see that 
x(S) = 2'~^x{^{p))- Let p = 2^™ be a cone radius of Z{p) at infinity as in 
Lemma 8.10. Note that m is polynomially bounded in the input size of S (given 
by the sparse bit size of the family of polynomials describing S). Consider the 
semialgebraic set T C -^n+r+m+i defined by 

p(x, y) = 0, = 2, - = 0, . . . , - = 0, \\xf + < z^. 

Clearly, T is homeomorphic to Z{p)p = Z{p) n + \\y\\'^ < p^}. Therefore, since 
p is a cone radius, we have x{Z{p)) = x{Z{p)p) = x*{Z{T)). By Theorem 7.1 we 
can compute x*{Z{T)) in FF^^. This implies that x{S) may be computed within 
the same resources. □ 



Remark 8.11 Theorem 8.9 easily extends to the case where we also allow inequal- 
ities h(x) > in the definition of the basic semialgebraic set. For instance, for 
5 = {x G R'' I p{x) = 0, h{x) > 0} consider 

Z := {{x, y) G R"+i | p{x) = 0, h{x) - y^}. 

The sets Zj^ := ^ H {y > 0} and Z_ := Z R {y < 0} are closed semialgebraic sets 
both homeomorphic to S and Z = Z+UZ-_. The formula x(Z+UZ_)-|-x(Z-|_nZ_) = 
x{Z+) + xi^-) then allows to compute x(>S') from the Euler characteristic of real 
algebraic varieties. 

8.3 Connected components and Betti numbers 

We are going to study here the following problems: 

#CC]g {Counting connected components) Given a semialgebraic set S, compute 
the number of its connected components. 

BETTi(fe)]R (kth Betti number of a real algebraic set) Given a real multivariate 
polynomial, compute the kth Betti number of its real zero set. 

BM-BETTi(fc)]R (fcth Borel-Moore Betti number of a real algebraic set) Given a 
real multivariate polynomial, compute the A;th Borel-Moore Betti number of 
its real zero set. 
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For the problems related to Betti numbers, we restrict the input to be a real 
algebraic set. Since we will only prove lower bounds for these problems, this re- 
striction makes our results stronger. Note that Betti(0)]r is just the restriction of 
#CCk to real algebraic sets. 

We will focus here on the discretized versions of the above problems, where the 
input polynomials have integer coefficients, and study these problems in the Turing 
model. 

The following upper bound was first shown by Canny [18] . 

Theorem 8.12 The problem #CC^ is in FPSPACE. 

Prom a result by Reif [60, 61] on the PSPACE-hardness of a generalized movers 
problem in robotics, it follows easily that the problem ^^CC^ is in fact FPSPACE- 
complete. We will give an alternative proof of the F PSPACE-hardness of this problem 
following the lines of [17]. This will also allow us to sharpen the lower bound 
by showing that #CC^ remains FPSPACE-hard when restricted to compact real 
algebraic sets. Based on this, we will prove the F PSPACE-hardness of the problems 
BETTl(fc)K and BM-Betti(/c)k. 

The following lemma follows by inspecting the usual NPR-completeness proof of 
Fease [9]. 

Lemma 8.13 For A £ there is a polynomial time Turing machine computing 
on input n G N a quantifier free first order formula G in the free variables 
xi, . . . , Xp(„) such that the projection 

{x G M^'^") I ^n{x) holds }^An M", (xi, . . . , ^ (xi, . . . , x„) 

is a homeomorphism. The inverse image of an integer point x G AD Z" is again 
integer and can be computed in polynomial time. 

Lemma 8.14 There is a polynomial time Turing machine computing from a quan- 
tifier free formula $ G J^jg in the free variables Xi , . . . , a polynomial /$ in 
Z[Xi,...,Xm,Yi,..., Fy(„)] such that the projection vr : -^R"',{x,y) ^ x 

induces for all e G { — 1, 1}^^™) a homeomorphism 

2{U) n {em > 0, . . . , eqirn)yq{m) > 0} ^ {x G M"* I ^{x) holds }. 

Proof. As in the NP^-completeness proof of FeaSr [9] the machine M performs 
the following (see also [19]). For each atomic formula of $ containing an inequality 
choose a new variable Y and replace 

p{X) > by p{X) - y2 = 
p{X) > by p{X)Y'^ -1 = 0. 
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In the resulting formula iteratively eliminate the connectives as follows: replace 

MUVi = by = 0, and A-=i^« = by ELiP? = 0. 

We end up with a single polynomial equation /$ = 0, which is easily seen to satisfy 
the claim of the lemma. □ 

Consider the following auxiliary problem: 

ReacHu {Reachability) Given real polynomials /, 5, decide whether there ex- 
ist points p G Z^n{f,g) and q G Z^n{f,h) which lie in the same connected 
component of Z^n(f). 

Proposition 8.15 The problem Reach^ is PSPACE-Lard. 

Proof. Assume L € PSPACE. In the proof of [17, Proposition 5.9] the configuration 
graph of a symmetric Turing machine deciding membership of 'u; G {0, 1}" to L was 
embedded in a certain way in Euclidean space as a compact onc-dimcnsional semi- 
linear set 5"^. More specifically, a polynomial time computable function mapping 
w G {0,1}" to i'TfmU'n,Vn) was constructed, where ^„ is a constant free additive 
circuit describing membership to C W^^^' , c is a polynomial, and Uyj , Vyj G 
{0, 1}^(") such that w £ LiS Uyj and Vyj are connected in Su,. Note that, in particular, 
the set A := {{w,x) G {0,1}" x W^"^^ | n G N,a; G Sy,} is contained in P°^^j and 
hence in P^. 

We apply Lemma 8.13 to the set A. Let G be the formula in the 

free variables Xi, . . . corresponding to the input size n -|- c(n) and let /„ G 

. . . , Yi, . . . , yg(„)] be the integer polynomial corresponding to ac- 

cording to Lemma 8.14. We know that can be computed from n in polynomial 
time. For w G {0,1}" let Hwi^w S IP^"^^ be the inverse images of {w,Uyj), (w,Vw) 
under the projection homeomorphism 

r„ := {x G W(^^ \ $„(x) holds } AnW+<^\ (xi, . . . , ^ (xi, . . . , 

Note that fiuj and t'^ are connected in r„ iff and are connected in Sw, which 
is the case iS w E L. 

According to Lemma 8.14, for any e G {—1,1}^^"^, the projection h- > x 
induces a homeomorphism 

^(/n) n {eiyi > 0, . . . , eq(n)yq{n) > 0} > T^. 

This implies that there exist points {f^w,il)j {^wiVi') ^ ^ifn) that are connected in 
Z{fn) iff fiw and Vy, are connected in T„. Define the integer polynomials g^j := 
fn{l^w,y),hw := gn{i^w,Y). Then u; G L iff the instance fn,9w,hw of the problem 
Reach^ has a solution. Moreover, fn,9w, hyj can be computed in polynomial time 
from w. □ 
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Remark 8.16 The proof of Proposition 8.15 shows that Reach^ remains PSPACE- 
hard when restricted to one-dimensional compact real algebraic sets. 

Lemma 8.17 For a compact Z CW^ let '^{Z) C M"+^ be the one-point compacti- 
Gcation ofZxR. Then we have 6^+i(S(Z)) = be{Z) for all leN. (This is also true 
for Z = $ with the convention that S(0) is a one point space.) 

Proof. The suspension S{Z) of a nonempty topological space Z is defined as the 
space obtained from the cylinder Z x [0, 1] over Z by identifying the points in each 
of the sets Z x {0} and Z x {1} obtaining the points vq and vi. Essentially, this 
is a double cone with basis Z and vertices vo,vi. It is well known that the Betti 
numbers of S{Z) and Z are related as follows (cf. [32, 55]): 

Assume, without loss of generality, that Z is nonempty. Since Z is compact, 
the one-point compactification ^{Z) of Z x M is homeomorphic to the suspension of 
S{Z) where the two vertices vq and vi of the double cone have been identified. This 
space is homotopy equivalent to the space obtained from the suspension S{Z) by 
connecting the vertices vq and vi with a one-dimensional cell. This space, in turn, 
is homotopy equivalent to the space obtained from S{Z) by attaching a circle at 
a point. Since this amounts to attach to Z only a cell we conclude that 



6^+1 (5(Z)) if£>0 
bi{S{Z)) + l if£ = 0. 



Combining this with (17), the claim fe^+i(E(Z)) = bi{Z) follows, for any £ G N. □ 

The one point compactification of a non-compact real algebraic set can be real- 
ized as a real algebraic set by a simple construction [10, p. 68]. For ^ G M" consider 
the homeomorphism (inversion with respect to the unit sphere with center ^) 
defined by 

Let / be a real polynomial of degree d with zero set Z C and assume that ^ ^ Z. 
Consider the polynomial := \\X - Cf^f{^ + \\X - ^\\~^{X - ^)) with zero set 
Z^ C M". If Z is unbounded then Z^ = t^{Z) U {^} is homeomorphic to the one- 
point compactification of Z. Note that if Z is empty, then Z^ consists just of the 
point 4- 

Theorem 8.18 For any /c G N both problems Betti(A;)^ and BM-Betti(A;)^ are 
FPSPACE-iiard with respect to Turing reductions. 
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Proof. Note first that the Borel-Moore and the usual Betti numbers coincide for 

compact sets. We denote by CBETTlg(A;) and CReach^ the restrictions of the 
problems Betti(A;)}^ and Reach^ to compact real algebraic sets. We know by 
Proposition 8.15 and Remark 8.16 that CReach^ is FPSPACE-hard. To prove the 
theorem, it is thus sufficient to establish a Turing reduction from CReach^ to 
CBetti^(A;). Our proof is similar to the one of [17, Lemma 5.20]. 

We first describe a Turing reduction from CBETTl}g(0) to CBETTl}g(fc), for fixed 
k > 0. Let the compact Z = Z{f) C M" be given by / G Z[Xi, . . . , Set 
/o := f + Co := (0, . . . , 0, 1) G M"+i and note that Co ^ ^(/o) = Z{f) x {0}. 

We recursively compute the sequence of polynomials /i, • • • as follows. Let 
1 < i < k and assume that G . . . , X^+j] has already been computed such 

that Ci_i := (0, . . . , 0, 1, . . . , 1) G 1^^^+* (n zeros, i ones) is not contained in 
Let fi-i denote the polynomial /i_ijnterpreted as a polynomial in Xi, . . . , 
where is a new variable and Ci_i := (Ci_i,0) G R"'+*+^. Note that = 

Z{fi-i) >^ We define now the polynomial fi := which results from 

fi-i by transformation with the inversion ij- w.r.t. the unit sphere with center 

si— 1 

Ci_i (sec the comments before Theorem 8.18). Note that = ^(Cj) ^ ^ifi) 

since - = 1 and Cj_i Z(/j_i). Then we have Z{fi) = S(2:(/j_i)) and 

Lemma 8.17 implies that bo{Z) = hk{Z{fk)). This gives the desired reduction from 
CBetti^(O) to CBetti^(A;). 

In order to show that CReach^ reduces to CBETTl{g(0) we first discuss an 
auxiliary construction. Assume we are given real polynomials f,g such that Z{f) C 
M" is compact and Z{f,g) is nonempty. Consider the one-point compactification 
Zf-g C 1^"+^ of the space Z{f) U {Z{f,g) x M). Topologically, this space is obtained 
from Z{f) by attaching a double cone with base Z{f,g) and identifying the two 
vertices of this cone. What is important is that all the points of Z{f, g) are connected 
in the new space. This is illustrated in Figure 1 below where Z{f) is the three closed 
curves, Z{g) is the dotted curve and, consequently, Z{f,g) is the four intersecting 
points. 
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Figure 1: An auxiliary construction. 

Using inversions as above, an equation of an algebraic set homcomorphic to Zj.^ can 
be easily computed from /, g. Let /i by a further polynomial such that ^(/, K) ^ 0. 
By attaching a double cone with basis Z{f, h) to we get a real algebraic variety 
^f;g,hi where all the points of Z{f,g) and Z{f,h), respectively, are connected. 

We describe now the Turing reduction from CReach^ to CBetti}^(0). For a 
given instance f,g,h& . . . , X„] of CReach^ we first check whether Z(f, g) or 

Z{f, h) is empty by two oracle calls. If this is the case, the corresponding reachability 
problem has no solution. Otherwise, we know that both Z{f,g) and Z{f,h) are 
nonempty. We compute now equations for the spaces Zf-g^h and Zf-gh (note that 
in the latter, all points of Z{g) U Z{h) have been connected). The spaces Zj-g^^ 
and Zf-gh have the same number of connected components iff there exist points 
p G Z{f,g) and q G Z{f,h) which lie in the same connected component of Z{f). 
Hence we get the desired reduction using two more oracle calls, one for Zf.g^h and 
one for Zf.gh- □ 

Remark 8.19 The Betti numbers modulo a prime p are defined similarly as the 
Betti numbers, but replacing the coefficient field Q by the finite field ¥p. It is 

easy to check that the proof of Theorem 8.18 also gives the FPSPACE-hardness of 
the computation of the A;th Betti number modp, and similarly for the Borel-Moore 
Betti numbers. 

9 Summary and final remarks 

We have summarized the results of this paper in Figure 2 which contains three 
diagrams showing results in the Turing model, over C, and over M. In this figure, an 
arrow denotes an inclusion, problems in square brackets are Turing-complete for the 
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class at their left, problems in curly brackets are many-one-complete for that class, 
and problems in angle brackets are hard for that class. The problems appearing in 
the figure are defined in the list below. Recall that if L denotes a problem defined 
over M or C, we denote its restriction to integer inputs by L". 

FPSPACE [#CC°] {BETTi(fe)S, BM-BETTi(fc)S> 



FpGCR [EULERg, EulerJ°] 

t 

ppGcc [Degree"] 

t 

FP#P 



GCR {#SASh" (CNF)} 
[#SASmO, #FeaSk°] 

GCC {#HNc°} 



FPARc 



FP^^'= [Degree] 



FPAR« 



FP** [Eulerk] 



#Pc {#HNc} #P. {#SAS, (CNF)} [#SASb,#Feask] 

Figure 2: Survey of main results. 

#FeaSr {Real algebraic point counting) Given a real multivariate polynomial, count the 
number of its real roots, returning oo if this number is not finite. 

#SASr {Semialgebraic point counting) Given a semialgebraic set S, compute its cardi- 
nality if S is finite, and return oo otherwise. 

#SASr (CNF) [Semialgebraic point counting) Given a semialgebraic set S in conjunctive 
normal form, compute its cardinality if S is finite, and return oo otherwise. 

Eulerr {Euler characteristic for basic semialgebraic sets) Given a basic semialgebraic 
set S decide whether S is empty and if not, compute x{^)- 

EuLERg {Modified Euler characteristic) Given a semialgebraic set S, decide whether it is 
empty and if not, compute its modified Euler characteristic. 

#CCr {Counting connected components) Given a semialgebraic set S, compute the num- 
ber of its connected components. 
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Betti(/;;)]k (fcth Betti number of a real algebraic set) Given a real multivariate polynomial, 
compute the fcth Betti number of its real zero set. 

BM-BETTl(fc)K (kth Borel-Moore Betti number of a real algebraic set) Given a real mul- 
tivariate polynomial, compute the fcth Borel-Moore Betti number of its real zero set. 

#HNc (Algebraic point counting) Given a finite set of complex multivariate polynomials, 
count the number of complex common zeros, returning oo if this number is not finite. 

Degree {Geometric degree) Given a finite set of complex multivariate polynomials, com- 
pute the geometric degree of its affine zero set. 

Other problems which appeared in this paper are listed below. The first three 
are NP]R-complete, the other two, NPc-complete. 

FeaSr {Polynomial feasibility) Given a real multivariate polynomial, decide whether it 
has a real root. 

SASk {Semialgebraic satisfiability) Given a semialgebraic set S, decide whether it is 
nonempty. 

DiMr {Semialgebraic dimension) Given a semialgebraic set S and d €N, decide whether 
dim 5 > d. 

HNc {Hubert's Nullstellensatz) Given a finite set of complex multivariate polynomials, 
decide whether these polynomials have a common complex zero. 

DiMc {Algebraic dimension) Given a finite set of complex multivariate polynomials with 
afiine zero set Z and d G N, decide whether dim Z > d. 

Remark 9.1 (i) To fix ideas, we assumed in the definition of the above problems 

that the input polynomials are given in sparse representation. However, note 
that choosing the dense encoding leads to polynomial time equivalent prob- 
lems. In order to see this, one just has to introduce additional variables that 
help to represent monomials of high degree by "repeated squaring" . The so- 
lution set of the new system of polynomial (in)cqualitics is homcomorphic to 
the original one. A similar remark applies for the encoding of polynomials by 
division free straight-line programs. 

(ii) Instead of restricting inputs to integer polynomials, one could allow also alge- 
braic (or real algebraic) coefficients with their standard binary encoding. The 
results in this paper would then hold as well and our proofs would only need 
some extra algorithmics, common in symbolic computation. 
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10 Open problems 



We believe that the developments in this paper open up a variety of meaningful new 
questions. To finish this paper we list some of them. 

Problem 1 Can one decide FeaSk in polynomial time with a black box for the 
Euler characteristic? 

Problem 2 It is known that the problem to count the number of connected com- 
ponents of a semialgebraic set is in FPAR]r. Is it hard in this class? We know that 
the corresponding result is true in the additive setting [17]. 

Problem 3 What is the complexity to check irreducibility of algebraic varieties 
over C? And what is the complexity of counting the number of irreducible compo- 
nents of algebraic varieties? 

Problem 4 Can Betti numbers of semialgebraic sets be computed in FPARr? We 
know that, in the additive setting, the computation of Betti numbers of semi-linear 
sets is FPARadd-completc [17]. 

Problem 5 What is the complexity to compute the multiplicity mu\tx{Z) of a 
point X in an algebraic variety Z7 And how about the computation of intersection 
multiplicities i{Z,A;x)'? 

Problem 6 What are the Boolean parts GCR and GCC of #P]^ and #Pc, respec- 
tively? 

Problem 7 Toda's theorem [68] states that PH C FP'^^. Is there an analogue of 
this over M or over C? 
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